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This  report  is  concerned  with  the  application  of  speckle  imagery 
techniques  to  the  task  of  space  object  identification  for  deep  space  objects. 

It  is  easy  to  see  that  because  of  the  very  large  range  involved,  a ground- 
based  imaging  system  would  have  to  have  a diameter  of  the  order  of  100  m 
to  provide  really  useful  resolution.  The  almost  certainly  impossible 
nature  of  a conventional  adaptive  optics  system  of  this  size  forces  us  to 
seek  other  approaches  to  the  imaging  problem.  The  approach  we  investi- 
gate in  this  report  is  based  on  a Michelson  interferometer  version  of  a 
speckle  imaging  system. 

The  technical  content  of  this  report  is  set  forth  in  two  chapters. 

The  first  of  these,  Chapter  2,  provides  a moderately  detailed  description 
of  the  Michelson  Interferometer  version  of  the  speckle  imaging  system. 

The  system  described  consists  of  three  units:  a stationery  4 m diameter 
telescope,  a mobile  1 m diameter  telescope,  and  a mobile  speckle -forming 
unit,  along  with  a flat  concrete  "apron1',  semi-circular  in  shape  and  with  a 
100  m radius.  During  an  extended  "exposure  period"  of  tens  of  minutes  to 
several  hours,  the  mobile  telescope  moves  about  the  apron,  and  together 
with  the  stationary  telescope,  provides  optical  inputs  to  the  speckle -forming 
unit,  based  on  which  a series  of  speckle  images  are  formed  and  recorded. 
After  the  "exposure  period,"  the  set  of  recorded  speckle  images  are  pro- 
cessed and  a high  resolution  image  of  the  object  in  question  is  developed. 
Chapter  2 describes  the  critical  features  of  each  unit,  and  the  system's 
mode  of  operation. 

In  Chapter  3,  the  problem  of  noise  in  speckle  imagery  and  speckle 
interferometry  is  addressed.  A rigorous  analytic  basis  is  developed  and 
formulas  are  obtained  for  the  rms  error  in  speckle  imagery  as  a function 
of  such  things  as  target  size,  target  brightness,  and  number  of  exposures. 
Sample  results  particularized  to  a full  circular  aperture  are  presented. 


These  will  be  generalized  to  a Michel  son  interferometer  type  aperture 
at  a later  time.  For  the  case  of  the  circular  aperture,  in  Eq.  (3. 161)  a 
simple  expression  is  given  for  the  required  target  brightness  so  that 
speckle  imagery  results  will  not  be  limited  by  photon  shot  noise,  but 
rather  by  the  limited  number  of  short  exposure  samples. 
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Introduction 


2.  1 

Space  object  identification  (SOI)  of  deep  space  objects  by  formation 
of  a visible  light  image,  from  a ground  station,  poses  a number  of  major 
problems.  The  two  most  demanding  considerations  derive  from  the  pro- 
blem of  imaging  through  atmospheric  turbulence  and  from  simple,  \/D  - 
type  diffraction  considerations.  For  deep  space  SOI  we  are,  of  necessity, 
forced  to  consider  aperture  diameters  of  the  order  of  10  m to  100  m.  Such 
aperture  sizes  are  almost  certainly  too  large  to  be  realized  with  conven- 
tional, full  aperture  telescopes,  and  also  so  large  as  to  probably  preclude 
consideration  of  adaptive  optics  in  the  manner  of  the  CIS-sensor,  as  the 
way  of  getting  around  atmospheric  turbulence  limitations.  We  believe, 
however,  that  it  is  possible  to  develop  the  desired  SOI  imaging  capability 
within  these  constraints,  and  shall  sketch  our  approach  in  the  following 
pages. 

The  key  to  the  achievement  of  the  desired  capability  lies  in  two 
techniques.  They  are,  first  the  exploitation  of  recent  developments  in 
speckle  imagery,  and  second  the  use  of  a continuously  adjustable  baseline 
Michelson  stellar  interferometer  type  telescope.  With  a reasonable  length 
data  collection  period,  appropriate  servo  systems  to  control  alignment  and 
path  length,  and  sufficient  data  processing  capability  this  type  of  system 
should  be  able  to  produce  high  resolution  visible  light  images  of  deep  space 
objects.  Because  of  the  low  apparent  angular  rate  of  deep  space  objects 
in  traveling  across  the  sky,  several  hours  should  generally  be  available  for 
recording  the  signals  necessary  to  allow  the  image  to  be  formed.  The 
extent  of  the  required  data  collection  period  need  not  be  a fundamental 
problem. 

In  this  chapter,  we  shall  be  concerned  with  the  development  of  an 
overview  description  of  this  type  of  deep  space  satellite  speckle  imager. 

We  shall  point  out  the  areas  where  we  foresee  important  technical  questions 


requiring  analysis,  and  will  indicate  the  critical  component  performance 
requirements.  Finally,  we  shall  indicate  the  type  of  performance  we  may 
expect  from  the  system. 


2.  2 


Basic  Considerations 


For  space  object  identification  by  inspection  of  imagery,  we  may 
take  as  our  basic  requirement  a resolution  of  6 = 15  cm  (s  6 inches). 

This  is,  of  course,  a somewhat  arbitrary  choice  — but  represents  what 
we  believe  is  necessary  if  we  are  to  be  able  to  do  more  then  simply  guess 
about  the  objects  identity  and  function  based  on  gross  features.  If  we 
identify  this  resolution  scale,  6 , with  the  spatial  frequency  half  wave- 
length at  half  the  apertures  cut-off  spatial  frequency,  £/\  , we  find  that 
for  an  object  at  range,  R , we  require 

[(i^A)'1  ]R  = 26 


or 


R \/J» 


U) 


Here  \ denotes  the  operating  wavelength,  which  we  may  take  to  be 
X = 0.  55  X lO^m,  while  £ denotes  the  maximum  separation  between 
points  on  the  collection  aperture.  For  a full  aperture  system  (i.  e.  , a 
simple  circular  aperture)  3 represents  the  aperture  diameter,  while 
for  a Michelson  stellar  interferometer  type  telescope,  3 denotes  the  max- 


For  imaging  of  an  object  in  a synchronous  orbit,  where  R es  3,6 
X 107  m,  the  required  aperture  extent  is 

b = RX/6 

= 132  m . (2) 

This  is  so  large  that  it  is  immediately  obvious  that  we  cannot  consider 
a full  aperture  system. 

Even  for  an  object  at  one-tenth  that  range  (i.  e.  , R = 3.  6 X 10^  m) 
the  required  aperture  diameter  (i.e.  , b - 13.2  m)  is  so  large  that  con- 
sideration of  a full  aperture  system  is  at  least  questionable.  Thus,  we 
may  conclude  that  we  must,  of  necessity,  contemplate  a synthetic  aperture 
type  system.  We  postulate  that  the  realization  of  this  synthetic  aperture 
system  will  be  based  on  the  Michelson  stellar  interferometer  type  telescope. 
In  the  balance  of  this  work  we  will  take  this  conclusion  for  granted.  The 
nominal  form  of  the  Michelson  stellar  interferometer  telescope  is  shown 
in  Figure  1 . 


figure  1.  Nominal  Form  of  Michelson  Stellar  Interferometer. 


The  sub  aperture  in  D and  the  baseline  site  is  b . 


Atmospheric  turbulence  will,  in  general,  limit  the  angular  resolution 
of  conventional  imaging  techniques  (whether  formed  with  a full  aperture 
telescope  or  a Michelson  stellar  interferometer  type  telescope)  to  6 5 

X 10*“®  rad,  or  larger.  This  resolution  is,  of  course,  incompatible  with 
the  desired  performance,  corresponding  to  a resolution  length  6 = 180  m 
at  a synchronous  orbit  range  of  R = 3.  6 X 107  m . The  use  of  adaptive 
optics  to  compensate  for  turbulence  effects  is  precluded  first  of  all  by  the 
very  weak  signal  strength  to  be  expected  in  the  self-referenced  adaptive 
optics  compensation  mode*,  but  perhaps  more  significantly,  by  the  fact 
that  available  wavefront  distortion  sensing  techniques  are  all,  in  essence, 
shearing  interferometer  or  Hartman  sensor  methods.  These  methods  are 
only  applicable  if  there  is  no  significant  gap  between  any  two  portions  of 
the  aperture.  But,  of  course,  for  the  Michelson  stellar  interferometer 
type  of  telescope  there  are  very  substantial  gaps  — so  that  we  have  no 
available  wavefront  distortion  sensing  technique  around  which  we  can  de- 
sign an  adaptive  optics  compensated  imaging  system. 

In  order  to  avoid  this  basic  limitation,  and  yet  circumvent  the 
atmospheric  turbulence  imposed  limitation  on  resolution,  we  plan  to  make 
use  of  speckle  imagery  techniques.  The  basic  speckle  interferometry 
technique  was  developed  by  Labeyrie1  , but  could  only  generate  a high 
resolution  correlogram  rather  than  the  target's  image  per  se.  The  cor- 
relogram  is  useful  for  interpretation  of  simple  objects  like  binary  stars 
and  resolved  stellar  disks,  but  is  totally  unsuited  to  our  objective  of  space 


* The  weakness  of  the  signal  strength  is  a direct  consequence  of  the  fact 
that  the  size  of  the  object  being  viewed  will,  in  general,  be  much  less 
than  the  turbulence  limited  resolution  size.  This  results  in  a very  low 
value  of  target  brightness  measured  in  stellar  magnitudes  per  resolution 
element,  and  it  is  this  quantity  that  determines  the  noise  limitation  in  the 
performance  of  self  referenced  adaptive  optics. 
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object  identification.  More  recently,  Knox  and  Thompson3  have  extended 
Labeyrie's  technique  to  allow  for  true  image  recovery  from  speckle  data 
of  even  very  complex  objects.  The  validity  of  this  technique  has  since  been 
demonstrated  in  computer  simulation  by  Nisenson,  Ehn,  and  Stacknik3. 

Thus  there  is  nothing  new  about  the  idea  of  using  speckle  imagery  techniques 
for  SOI.  Our  particular  inovation  lies  in  the  concept  of  applying  speckle 
imagery  techniques  to  data  taken  with  a variable  baseline  Michelson  stellar 
interferometer  type  telescope,  and  in  the  particular  form  of  this  optical 
configuration. 

Because  of  our  plans  to  use  a Michelson  stellar  interferometer,  a 
very  large  number  of  short  exposures  will  have  to  be  formed.  Of  the  order 
of  10s  to  103  will  be  required  for  each  baseline  setting.  This  is  needed 
to  provide  adequate  statistical  averaging.  Moreover,  since  the  interfero- 
meter can  handle  only  one  limited  set  of  spatial  frequencies  at  a time,  a 
large  number  of  interferometer  settings  will  be  required.  We  shall  suggest 
a design  for  which  about  1000  settings  are  required  for  observation  of  a 
satellite  in  a synchronous  orbit.  Thus  between  10^  and  10^  total  short 
exposures  will  have  to  be  generated  with  each  exposure  lasting  about  10 "2 
sec,  this  will  require  between  about  20  minutes*  and  3 hours  for  the  total 
image  recording  process.  For  a satellite  in  a lower  orbit  because  a smaller 
maximum  value  of  the  baseline,  £ , is  required  to  provide  the  desired  reso- 
lution, a lesser  number  of  interferometer  baseline  setting  is  needed,  and  so 

the  total  image  recording  time  will  be  smaller.  The  time  should  vary  about 

« 

as  the  square  of  the  satellite  altitude  or  range. 

Because  we  intend  to  rapidly  vary  the  interferometer  baseline  it  will 
be  necessary  to  introduce  servo  controls  to  handle  alignment  and  equalization 
of  the  path  length  of  the  two  arms  of  the  interferometer.  The  path  length 


* This  ignores  time  required  to  vary  the  interferometer  baseline  setting. 


equalization  servo  will,  of  necessity,  be  based  on  white  light  fringe  sensing. 
If  this  servo  is  run  with  a high  enough  servo  bandwidth,  it  will  not  only  pro- 
vide basic  alignment  control,  but  will  also  reduce  the  turbulence  imposed 
limitation  on  the  allowable  spectral  bandwidth.  This  should  have  a major 
impact  on  the  required  short  exposure  time. 

With  this  introduction  and  consideration  of  basic  matters  we  are 
now  ready  to  take  up  the  matter  of  providing  a description  and  definition 
of  the  interferometer  speckle  imaging  system.  This  is  treated  in  the 
next  section. 

2.  3 System  Description 

The  system  consists  of  four  basic  components.  The  first  of  these 
is  a very  flat  semicircular  concrete  surface  with  a radius  L = 132  m, 
which  we  shall  refer  to  as  the  "apron.  " The  second  component,  located 
at  the  center  of  the  circle,  of  which  the  semicircle  is  a part,  is  a D = 4 m 
diameter  telescope,  which  we  shall  refer  to  as  the  "base  aperture."  The 
third  component  of  the  system  is  a d = 1 m diameter  telescope  mounted  on 
a mobile  base  and  able  to  position  itself  anywhere  on  the  apron.  We  refer 
to  this  as  the  "mobile  aperture.  " The  fourth  and  last  of  the  basic  compon- 
ents is  a mobile  optical  system  to  which  both  the  base  aperture  and  the 
mobile  aperture  relay  the  radiation  received  from  their  field-of-view.  In 
this  unit,  which  we  call  the  "speckle  former,  " the  beams  from  the  two 
apertures  are  combined  to  form  the  desired  speckle  pattern. 

The  speckle  former  unit  contains  the  basic  servo  capabilities  that 
result  in  the  images  formed  by  the  two  apertures  being  properly  registered 
with  respect  to  each  other.  It  also  contains  the  basic  white  light  fringe  path 
length  equalization  servo.  The  two  aperture  units,  i.  c.  , the  base  aperture 
unit  and  the  mobile  aperture  unit  each  have  pointing  servos  to  allow  them  to 
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track  the  target  being  observed  and  to  direct  their  output  beams  to  the 
entrance  ports  of  the  speckle  former  unit.  These,  however,  are  rela- 
tively modest  performance  servos  compared  to  those  in  the  speckle  former 
unit.  It  is,  in  fact,  the  angular  alignment  (i.  e.  , the  two  beam  registration) 
and  the  path  length  equalization  Bervos  in  the  speckle  former  unit  that  will 
produce  the  necessary  degree  of  optical  alignment  tolerance  needed  to  allow 
the  two  arms  of  the  interferometer  to  work  together,  in  a quasi  coherent 
manner*.  Because  the  mobile  aperture  unit  is  continuously  moving,  the 
performance  of  these  servos  is  critical  to  the  ability  of  the  system  to 
function. 

The  angle  tracking  alignment  achieved  in  the  speckle  forming  unit 
is  merely  good  enough  to  keep  two  low  resolution  (turbulence  limited)  spots 
reasonably  well  superimposed  on  each  other.  The  path  length  equalization 
part  of  the  speckle  forming  units  servo  control  is  required  to  insure  that 
the  optical  path  lengths  from  the  source  to  the  speckle  forming  unit's  focal 
plane  via  the  base  aperture  and  via  the  mobile  aperture  are  equal.  To 
accomplish  this  the  speckle  forming  unit  has  to  be  positioned  part  way  bet- 
ween the  base  aperture  and  the  mobile  aperture.  (Like  the  mobile  aperture, 
the  speckle  forming  unit  will  have  to  be  continuously  moving  on  the  apron.  ) 
If  the  target  were  directly  overhead,  the  speckle  forming  unit  would  have  to 
be  located  exactly  half  way  between  the  two  apertures,  but  with  a target 
located  away  from  the  zenith  the  speckle  forming  unit  would  have  to  be 
positioned  closer  to  one  aperture  than  to  the  other.  In  all  cases  the 
criteria  for  the  speckle  forming  unit  to  be  in  a satisfactory  position  is  that 
white  light  interference  fringes  can  be  formed  in  its  focal  plane.  Obviously, 


* Actually,  because  of  the  nature  of  speckle  imagery,  the  quasi  coherence 
requirement  is  considerably  less  severe  than  we  would  have  if  we  required 
full  coherence. 


it  will  not  be  possible  to  control  the  speckle  forming  units  position  with 
sufficient  accuracy  to  insure  this.  Accordingly,  a fine  path  length  ad- 
justing element  will  be  incorporated  into  the  optical  train  within  the 
speckle  forming  unit.  This  path  length  adjustment  element  will  be  servo 
controlled  using,  as  a basis  for  generation  of  a servo  error  signal,  a type 
of  white  light  fringe  visibility  sensor.  This  will  not  only  compensate  for 
geometric  path  length  effects,  but  will  also  accommodate  the  grosser 
aspects  of  the  turbulence  induced  path  length  variations.  This  will  sub- 
stantially increase  the  allowed  spectral  bandwidth  — from  a nominal  value 
of  il  ^ (r0/L)X.  ^(0,  1/132)  5000  X = 3. 8 A based  on  the  turbulence  effects 

over  the  full  interferometer  baseline  to  AX  (r0/D)X  ^ (0.  1/4)  5000  A = 
o 

125  A , which  value  is  set  by  the  turbulence  effects  across  the  larger  sub 
aperture  only.  The  key  to  this  improvement  is  that  the  path  length  servo 
in  the  speckle  forming  unit  be  fast  enough  to  track  the  large  scale  (and 
therefore,  relatively  slowly  changing)  turbulence  variations. 

To  facilitate  visualization  of  the  overall  system,  in  Figure  2 we 
have  depicted  the  several  units  making  up  the  system,  as  they  might  ap- 
pear on  the  apron.  In  the  following  sections  we  describe  some  of  the 
more  significant  features  of  each  unit.  We  will  then  describe  the  more 
subtle  aspects  of  the  system's  principles  of  operation,  and  finally  will 
provide  some  estimates  of  the  expected  performance. 
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Figure  2.  Overview  uf  J>ynlh«\ic  Aperiute  Sp»tWU  Imagery  SyMvtu. 

l'lie  «pron  »•  a very  flat  concrete  mrUce  covering  « »«nucircU  of 
Ii2  m radiue.  The  mobile  aperture  unit  h*i  an  aperture  of  1 m 
diameter  while  the  b*»e  aperture  unit  has  an  aperture  of  4 m diameter. 

The  speckle  unit  is  so  positioned  that  the  optical  path  length  from  the  source 
to  the  speckle  unit  via  the  two  collection  apertures  are  essentially  equal. 


2.4  Base  Aperture  Unit 

The  base  aperture  unit  consists  of  a D = 4 m diameter  telescope 
(in  a fixed  position  dome),  with  auxilary  optical  and  servo  units.  The 
exact  form  of  the  telescope  mount  (i.  e.  , altitude  inside  azimuth,  polar 
axis,  etc.)  is  relatively  immaterial.  What  is  important  is  that  the  tele- 
scope aperture  have  only  a rather  small  central  obscuration  (much  less 
than  one  meter  in  diameter),  and  that  by  some  means  (nominally  a hollow- 
axis  gimbal  system)  the  telescope  focus  be  brought  off  the  mount  so  its 
position  is  stationary,  i.  e.  , independent  of  the  telescope's  orientation. 

The  converging  beam,  just  before  it  reaches  the  focal  plane  is  beam  split 
by  a dichroic  filter,  splitting  off  about  a 1000^  portion  around  \ = 5500  A 
for  use  in  forming  the  speckle  image  and  in  generating  control  signals  for 
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the  speckle  image  forming  unit.  The  portion  of  the  beam  passed  by  the 
dichroic  beam  splitter  is  recollimated  to  a diameter  C is  0.2  m . (The 
recollimating  optics  are  refractive  so  no  additional  aperture  central 
obscuration  is  introduced  at  this  point. ) 

Because  a number  of  reflections  are  involved  in  bringing  the  beam 
off  the  telescope  mount,  the  beam  is  subject  to  rotation  about  its  own  axis, 
the  exact  amount  of  rotation  being  dependent  on  orientation  of  the  telescope. 
To  remove  this  rotation  the  recollimated  beam  (D  = 0.2  m diameter)  is 
passed  through  a K-mirror  assembly  which  is  driven  open  loop,  an  amount 
calculated  from  the  telescope  orientation.  The  accuracy  of  the  K-mirror 
setting  need  only  be  good  enough  to  remove  rotation  effects  to  about  ± 1°  . 
(Fine  rotation  corrections  will  be  made  later  on  a servo  controlled  K- 
mirror  device  in  the  speckle  forming  unit. ) 

Finally,  the  recollimated  beam  is  reflected  off  of  a two  axis  beam 
steering  unit,  and  transmitted  to  a beam  receiver  aperture  of  the  speckle 
forming  unit.  To  insure  accurate  positioning  of  the  beam  on  the  receiver 
aperture  of  the  speckle  forming  unit,  an  infrared  "pilot  beam"  is  trans- 
mitted from  that  receiver  aperture  and  sensed  by  a unit  that  looks  through 
the  two  axis  beam  steerer,  and  also  views  the  direction  of  the  recollimated 
beam. 

The  general  layout  of  the  base  aperture  unit  is  indicated  in  Figure  3. 


- 14  - 


!«|  infrared  pilot  beam 
'Jr  ^((rom  entrance  port 
{ { of  speckle  forming  unit) 

Figure  3.  Optical  Block  Diagram  for  Baae  Aperture  Unit. 

2.  5 Mobile  Aperture  Unit 

The  mobile  aperture  unit,  considered  as  an  optical  train  will  be 
essentially  the  same  as  the  base  aperture  unit,  except  for  the  fact  that 
the  telescope  aperture  diameter  will  be  d = 1.  0 m,  and  the  recollimated 
beam  diameter  will  be  d = 0.  05  m . This  beam  will  also  be  transmitted 
to  the  speckle  forming  unit,  but  will  enter  through  its  own  receiver  port. 
The  mobile  aperture  unit  will  have  the  same  sort  of  recollimation  optics, 
dichroic  beam  splitter,  telescope  pointing  control  sensor,  K-mirror,  and 
two  axis  beam  steering  mirror  as  the  base  aperture  unit.  Of  course,  they 
will  be  matched  to  a d = 0.05  m,  (rather  than  a D = 0.2  m)  recollimated 
beam  diameter. 

Aside  from  telescope  aperture  and  recollimated  beam  diameter 
differences,  the  most  substantial  differences  between  the  base  aperture 
unit  and  the  mobile  aperture  unit,  is  the  fact  that  the  mobile  aperture  unit 
will  be  mounted  on  some  sort  of  vehicle  and  will  be  able  to  move  about  on 


the  apron.  In  operation  we  plan  to  have  the  mobile  aperture  unit  con- 
tinuously moving,  at  a low  rate.  Nominally  the  rate  will  be  about  0.  1 
to  1.0  m/sec  . This  motion  and  associated  "instabilities"  will  impose 
servo  bandwidth  requirements  on  all  of  the  servo  loops  in  the  mobile 
aperture  unit,  but  at  this  point,  this  need  be  of  no  particular  consider- 
ation to  us  here. 

2.  6 Speckle  Forming  Unit 

The  central  and  most  critical  part  of  the  entire  system  is  the 
speckle  forming  unit.  It  is  here  that  the  beams  from  the  base  aperture 
unit  and  the  mobile  aperture  unit  are  brought  together  aligned  with  re- 
spect to  each  other  in  angle  and  in  path  length,  and  then  combined  on  a 
common  focal  plane  to  form  the  desired  image.  The  speckle  forming  unit, 
like  the  mobile  aperture  unit,  is  mounted  on  wheels  and  moves  on  the 
apron. 

Its  position  is  always  part  way  between  the  base  aperture  and  the 
mobile  aperture,  the  nominal  position  being  controlled  in  accordance  with 
the  geometric  considerations  indicated  in  Figure  4,  to  equalize  the  path 
length  via  the  two  apertures.  The  speckle  forming  unit  utilizes  a laser 
radar  unit  to  measure  the  distance  to  the  base  aperture  unit  and  to  the 
mobile  aperture  unit,  but  only  to  get  into  approximately  the  correct  position. 
Fine  adjustment  of  the  path  length  is  made  on  the  beam  from  the  mobile 
aperture  unit,  after  it  has  entered  the  port  in  the  speckle  forming  unit. 

The  adjustments  are  accomplished  with  a device  of  the  form  shown  in 
Figure  5,  which  we  refer  to  as  the  optical  path  expander. 

The  signal  to  control  the  optical  path  expander  setting  is  obtained 
from  a quasi  white  light  interferometer  unit.  The  unit  combines  equal 
diameter  portions  of  the  beams  from  the  base  aperture  unit  and  from  the 
mobile  aperture  unit  by  means  of  a beam  splitter.  The  diameter  of  each 
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portion  will  correspond  to  a diameter  of  about  ^ r0  , i.  e. , about  0.  05  m 

back  at  the  collection  apertures.  Before  being  combined  at  the  50:50  beam 

splitter /combiner,  one  path  will  be  phase  modulated  with  a piezoelectrically 

o 

driven  mirror,  about  ± 500  A , at  a modulation  frequency  of,  say,  3 x 10a 

Hz.  After  being  combined  the  composite  beam  will  be  reflected  off  of  a 

diffraction  grating  and  then  focused  on  an  intensified  linear  detector  array, 

o 

with  the  individual  detectors  having  a spectral  resolution  of  AX  = 1 A , 
and  there  being  a total  of  1000  detectors.  The  optical  configuration  is 
indicated  in  Figure  6. 
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Figure  4.  Speckle  Forming  Unit  Positioning  Consideration*. 

A whit*  light  (ring*  condition  null  b*  aatiafiad  in  ord*r  for  * uaeful  iptckl, 
pattern  to  b*  formed.  To  »chl*v*  thi*  th*  apeckl*  forming  unit  mint  b*  portioned 
,o  that  AB  + CD  m A'B'  ♦ C'D*  , who-  • th*  approximation  r*fl*ct*  th*  fact  that 
allowanc*  mint  b*  mad*  for  optical  pain  length*  within  th*  baa*  apertur*  unit,  th* 
mobile  aperture  unit,  and  th*  apeckl*  forming  unit. 


Fixed  Mirror  Pair 


Figure  5.  Optical  Path  Expander. 


The  optical  path  length  from  input  to  output  can  be 
varied,  (without  effecting  any  other  property  at  the  optical 
beam)  by  varytr.  the  separation  between  the  lixed  mirror 
pair  and  the  movable  mirror  pair.  It  is  perhaps  worth  noting 
that  the  output  beam  will  remain  parallel  to  the  input  despite 
a small  lilt  of  either  mirror  pair. 
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Figure  6.  White  Light  Interferometer  Optical  Layout. 
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In  operation  the  optical  path  expander  setting  will  be  slowly  moved 

through  a path  length  variation  range  of  perhaps  as  much  as  1.0  m,  though 

more  likely  a search  over  a range  of  0.  1 m should  be  sufficient.  For  the 
o 

= 1 A spectral  resolution  at  the  spectrally  dispersed  focal  plane,  when 
the  path  expander  setting  makes  the  paths  nearly  equal  (in  this  case  to 
within  about  ± 1.25  mm),  the  individual  detectors  will  see  intensity  oscil- 
lations at  the  3 kHz  phase  modulation  frequency.  The  exact  setting  of  the 
path  length  expander  can  then  be  refined  by  comparing  the  phase  of  the 
modulation  signal  as  seen  by  each  of  the  detectors.  With  suitable  data 
processing  the  path  lengths  can  eventually  be  equalized  to  within  a fraction 
of  a wavelength,  in  which  case,  all  of  the  detectors  will  see  inphase  mod- 
ulation of  3 kHz,  and  the  sum  of  the  detector  outputs,  i.  e.  , the  composite 
output  signal  will  show  strong  modulation  at  this  frequency. 

The  initial  alignment  of  the  optical  path  length  expander  will  be 

performed  with  both  the  mobile  aperture  unit  and  the  speckle  forming  unit 

stationary.  Once  the  white  light  fringe  condition  is  achieved  the  two  units 

will  be  allowed  to  move.  The  path  length  expander  device's  setting  will  be 

controlled  by  means  of  a relatively  slow  dither  of  the  setting,  say  at  200  Hz 
o 

over  ± 200  A . It  will  be  observed  that  the  depth  of  modulation  of  the  3 kHz 
variation  in  the  composite  signal  will  oscillate  at  a 200  Hz  rate.  The  nature 
of  this  oscillation  will  indicate  the  magnitude  and  direction  of  the  required 
correction  of  the  optical  path  expander's  setting,  thus  providing  a basis 
for  closed  loop  control  of  the  path  length  equalization.  To  the  extent  that 
the  optical  path  expander  setting  is  away  from  the  middle  of  its  range, 
this  can  be  used  to  fine  tune  the  motion  of  the  speckle  forming  unit  as  it 
follows  the  motion  of  the  mobile  aperture  unit. 

With  the  path  lengths  equalized  the  main  portion  of  the  speckle 
forming  unit  can  now  come  into  play.  To  form  the  speckle  pattern,  the 
optics  "simulates"  the  desired  portion  of  a L = 132  m diameter  aperture, 


f / 5 0 optical  system.  The  desired  portions  of  the  aperture  are  two 
circular  regions  corresponding  to  D = 4 m diameter  and  d = 1 m di- 
ameter with  a separation  corresponding  to  the  current  separation,  st  , 
of  the  base  aperture  unit  and  the  mobile  aperture  unit.  It  is  a simulation 
of  these  optics  in  the  sense  that  the  actual  entrance  aperture  is  only  2 = 

1.  32  m,  with  an  effective  focal  length  of  3 = 66  m.  The  two  beams  are 
first  reduced  from  their  initial  diameters  of  D = 0.2  m and  d = 0. 05  m 
to  D7  = 0.04  m and  d7  = 0.01  m.  They  are  then  put  into  two  periscope 
arms,  each  of  which  can  pivot  about  a point  at  the  center  of  their  injection 
apertures.  The  output  of  the  periscope  arms  are  above  a 2 = 1.  32  m di- 
ameter mirror,  which  with  its  secondary  optics,  has  an  effective  focal 
length  of  3 = 66  m.  The  output  of  the  two  periscopes  is  collected  by  this 
mirror  optics  and  brought  to  a focus.  If  the  two  spots  are  properly  ori- 
ented and  positioned  so  that  they  overlap*,  they  will  produce  a speckle 
pattern  matching  that  which  would  have  been  formed  by  the  corresponding 
portions  of  the  aperture  of  an  L = 132  m aperture  diameter,  f/50  telescope! 

In  order  to  insure  that  the  two  beams  are  properly  registered,  just 

before  the  two  beams  reach  the  focal  plane  they  will  pass  through  a dichroic 

filter  which  will  pass  the  narrow  band  protion  (about  200  X)  to  be  used  in 

forming  the  speckle  image  pattern,  but  will  direct  the  principal  part  of  the 

o 

energy  (i.  e.  , the  rest  of  the  1000  A bandwidth)  to  an  imager/tracker  section. 
Using  high  speed  chopping  to  block  first  one  and  then  the  other  of  the  two 
beams,  the  imager /tracker  will  be  able  to  sense  the  positioning  error  of 
each  of  the  two  spots  and  generate  a beam  direction  error  signal.  This 


* It  should  be  noted  that  the  individual  spot  sizes  will  be  of  the  order  of 
X/r0  <%s  0.  5 X 10"®/0.  1 = 5 X 10"®  rad.  So  long  as  the  registration  error 
between  the  two  spots  is  random  and  small  compared  to  X/r0  , it  can  be 
considered  as  equivalent  to  a small  increment  in  the  wavefront  tilt  portion 
of  the  optical  effects  of  turbulence,  and  will  be  accommodated  in  the  speckle 
image  data  processing. 


r 
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signal  can  be  generated  at  a high  data  rate  and  used  to  control  a two  axis 
fine  adjustment  pointing  units  at  the  input  to  each  of  the  periscopes.  The 
data  rate  should  be  sufficient  to  allow  a servo  bandwidth  that  can  track  out 
any  residual,  high  frequency  pointing  error  due  to  the  motion  of  the  mobile 
aperture  unit  or  of  the  speckle  forming  unit. 

In  addition,  in  order  to  correct  for  any  residual  rotation  about  the 
beam  axis  error,  the  imager/tracker  will  look  for  any  secondary  objects 
in  the  field-of-view  such  as  a star  which  is  momentarily  near  the  object, 
and  will  determine  the  K-mirror  rotation  correction  required  to  make  the 
two  beams  match.  This  will  be  implemented  by  a fine  control  K-mirror 
corrector  located  just  after  the  port  for  the  beam  from  the  base  aperture 
unit.  The  data  rate  for  control  of  this  K-mirror  will  probably  be  fairly 
low,  but  since  the  required  rotation  correction  should  only  change  rather 
slowly,  this  should  prove  adequate. 

In  Figure  7 we  indicate  the  optical  train  associated  with  the  speckle 
forming  unit  from  the  two  entrance  ports  to  the  two  periscopes  to  the  speckle 
and  the  tracker/imager  focal  planes. 

With  the  nature  of  the  basic  units  thus  defined  we  are  now  ready 
to  consider  the  operational  procedure  and  the  nominal  performance 
parameters.  We  take  this  up  in  the  following  sections. 

2.  7 Operational  Procedure 

When  a target  of  interest  reaches  a position  in  the  sky  where  we 
wish  to  start  recording  its  speckle  image  the  mobile  aperture  unit  will  be 
moved  to  within  four  meters  of  the  base  aperture  unit.  The  speckle  form- 
ing unit  will  be  positioned  near  the  two  aperture  units,  and  so  located  and 
operated,  that  it  can  receive  the  two  reduced  diameter  beams  sent  to  it 
from  the  base  aperture  unit  and  the  mobile  aperture  unit. 
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Figura  7.  Optical  layout  of  tha  .Speckle  Forming  Unit. 


At  this  point,  the  two  telescopes  will  be  pointed  toward  the  target 
and  start  to  track  it.  The  required  tracking  precision  will  be  5 prad, 

(2-tf)  . After  passing  the  collected  beam  in  each  aperture  unit  through  a 
field- stop  and  a dichroic  beam  splitter  (which  stripped  out  the  spectral 
portion  of  the  beam  to  be  used  for  tracking),  the  portion  of  the  beam  to  bo 
sent  to  the  speckle  forming  unit  will  be  recollimated  and  directed  to  the 
appropriate  entrance  port  of  the  speckle  forming  unit.  There  the  two 
beams  will  be  combined  at  a 50:50  beam  splitter,  after  one  of  the  beams 
has  gone  through  the  path  length  expander.  (A  selected  spectral  portion 
of  each  of  the  two  beams  will  have  first  been  split  off  to  be  used  in  form- 
ing the  speckle  image.)  The  combined  beams,  "joined"  at  the  50:50  beam 
splitter /combiner,  will  be  reflected  off  a diffraction  grating  and  sent  to  an 
array  of  detectors,  each  detector  corresponding  to  a narrow  spectral  band. 
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The  path  length  expander  will  be  slowly  scanned  through  its  control 
range.  Because  one  of  the  two  beams  going  into  the  50:50  beam  splitter/ 
combiner  was  subject  to  a 3 kHz  phase  modulation,  when  the  two  path 
lengths  are  nearly  equal,  i.  e.  , the  path  length  expander  is  near  its  pro- 
per setting,  the  signal  on  each  detector  will  be  observed  to  have  a signi- 
ficant 3 kHz  modulation.  At  this  point,  using  the  small  amplitude,  200  Hz 
modulation  of  the  path  length  expander,  a servo  loop  is  closed  to  make 
the  path  lengths  equal.  At  the  same  time  portions  of  the  two  beams  going 
to  the  speckle  imaging  plane  have  been  diverted  to  the  imager/tracker 
sensor  where  control  signals  were  generated  to  get  the  angular  alignment 
of  the  two  received  beams  adjusted,  so  that  on  the  speckle  image  plane 
the  two  are  registered  on  each  other  to  within  a small  fraction  of  \/r0  . 

At  this  point,  the  system  will  be  almost  ready  to  start  recording  speckle 
images. 

Before  the  system  is  fully  ready  to  start  the  recording  of  speckle 
images,  the  beam  rotation  will  have  to  be  adjusted  through  the  K-mirror 
assembly  in  the  speckle  image  unit.  This  will  be  done  by  slowly  nodding 
first  one  telescope  and  then  the  other  along  the  great  circle  that  passes 
through  the  zenith,  while  continuously  tracking  the  target.  This  will 
cause  the  image  as  seen  at  the  imager/tracker  to  similarly  nod.  (The 
nodding  amplitude  need  only  be  a few  arc  seconds.)  The  K-mirrors  will 
then  be  adjusted  so  that  both  beams  produce  images  that  nod  back  and 
forth  along  the  same  (or  parallel)  paths.  With  the  K-mirrors  thus  adjusted, 
the  system  will  be  fully  ready  to  record  short  exposure  speckle  images. 

The  individual  short  exposures  will  be  allowed  10  msec  each. 

During  any  one  exposure  a range  of  spatial  frequencies  is  recorded  cor- 
responding to  a 4 m diameter  range,  i.  e. , to  a range  of  |ah|  = D/X  = 

4/(.  5 x 10"®)  = 8 X 10s  . The  mobile  aperture  unit  will  be  continuously 
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moving  so  that  a slightly  different  set  of  spatial  frequencies  is  recorded 
in  each  frame.  Because  we  desire  between  10s  and  103  exposures  con- 
taining each  spatial  frequency,  the  speed  of  the  mobile  aperture  unit  will 
correspond  to  4 m displacement  in  10a  to  103,  10  msec  exposures,  i.  e.  , 

4 m/ sec  to  0.4  m/sec.  With  this  speed  the  mobile  aperture  unit  will 
traverse  a semicircle  about  the  base  aperture  unit,  and  then  increment 
its  distance  from  the  base  aperture  unit  by  4 m and  repeat  its  semicir- 
cular motion.  At  the  end  of  each  semicircular  traverse  the  distance  to 
the  base  aperture  unit  will  be  incremented,  until  finally  the  last  traverse 
will  be  made  at  a separation  between  the  mobile  aperture  unit  and  the  base 
aperture  unit  that  corresponds  to  the  desired  resolution.  This  will  depend 
principally  on  the  range  to  the  target,  the  total  traverse  time,  and  path 
length  being  proportional  to  the  range  to  the  target. 

In  practice  it  may  prove  convenient  to  move  the  mobile  aperture 
unit  at  twice  (or  three  times)  the  nominal  speed  and  increment  the  radius 
of  the  semicircle  by  one-half  (or  one-third)  of  the  nominal  4 m increment. 
This  will  reduce  the  number  of  short  exposures  containing  any  one  spatial 
frequency  during  a single  semicircular  traverse,  but  will  cause  that 
spatial  frequency  to  be  contained  in  speckle  images  recorded  on  other 
traverses. 

During  the  system's  operation  it  will,  of  course,  be  necessary  not 
only  to  record  each  speckle  image,  but  also  to  record  the  location  of  the 
mobile  aperture  unit  since  this  will  determine  the  actual  spatial  frequency 
content  of  the  image.  In  addition,  it  will  be  necessary  to  develop  the 
speckle  imagery/turbulence  transfer  function  for  each  spatial  frequency. 
ThiB  can  be  obtained  by  taking  advantage  of  the  periodic  presence  of  a star 
in  the  field-of-view  near  the  target  object.  At  14™  magnitude  there  are  an 
average  of  250  stars  per  square  degree.  Thus,  whenever  the  target's 
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angular  position  moves  by  one  degree  relative  to  the  celestial  background, 
approximately  2.87  stars  will  have  passed  within  ± 100  prad  of  the  target 
object.  The  speckle  forming  unit  will  be  designed  to  split  off  the  image  of 
such  a star  and  form  a speckle  image  on  a separate  focal  plane.  In  this 
way  speckle  calibration  data  can  be  obtained  for  various  spatial  frequencies*. 
There  will,  of  course,  be  sizable  gaps  in  the  spatial  frequency  calibration 
data  because  of  the  uncontrolled  nature  of  the  separation  between  the  base 
aperture  unit  and  the  mobile  aperture  unit  at  the  time  when  a useful  star 
moves  into  the  field-of-view,  but  because  of  the  smoothness  of  the  cali- 
bration function  versus  spatial  frequency,  there  should  be  no  problem  in 
interpolating  to  fill  in  the  missing  data. 

2.  8 Total  Exposure  Time 

As  noted  initially,  the  total  exposure  time  and  the  distance  the 
mobile  aperture  unit  will  have  to  travel  is  proportional  to  the  square  of 
the  targets  range.  For  a target  at  the  maximum  range  of  3.  6 X 107m  , 
we  can  estimate  the  mobile  aperture  travel  distance  as 

s = | tt 

where  Jfr  = 132  m is  the  maximum  separation  between  the  base  aperture 
unit  and  the  mobile  aperture  unit,  and  A = 4 m is  the  radial  distance  bet- 
ween the  semicircular  paths  traversed  by  the  mobile  aperture  unit.  Thus 
we  get 


s = |tt(132)»/4  = 6.842  x 10s  m . 


* Some  method  of  assuring  that  the  star  observed  is  not  actually  a binary 
star,  or  of  allowing  for  that  possibility,  will  have  to  be  implemented. 


Lj 
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With  a traverse  speed  of  v = 0.4  m/sec  (corresponding  to  1000  exposures 
of  each  spatial  frequency)  the  total  exposure  time  will  be 


T = s/v  = 6. 842  X 103/0. 4 

= 1. 711 X 104  sec 

= 4. 75  hrs 

We  must,  of  course,  add  to  this  some  time  for  the  change  from  one  semi- 
circle radius  to  the  next  semicircle  radius,  and  some  time  for  initial 
alignment. 

The  total  time  required  for  high  resolution  photography  of  a satel- 
lite in  a synchronous  equatorial  orbit  (range  equals  3.6  X 107m),  is  thus 
seen  to  be  rather  small.  If  the  range  to  the  satellite  were  smaller,  say 
5 X 10P  m (or  1 X 107  m)  the  required  exposure  time  would  drop  to  an  even 
smaller  value  from  4.  75  hours  to  only  5.  5 minutes  (or  22  minutes).  If  the 
100  exposures  rather  than  1000  exposures  per  spatial  frequency,  proves  to 
be  acceptable,  then  the  exposure  time  could  be  reduced  by  another  factor  of 
ten,  to  a few  minutes  or  even  fractions  of  minutes. 


2. 9 Conclusions 


Based  on  the  general  concept  definition  provided  here  and  the  very 
preliminary  analysis  which  we  have  gone  through,  we  conclude  that  it 
should  be  possible  to  use  speckle  techniques  to  form  high  resolution  images 
of  objects  in  deep  space.  The  components  required  are  reasonably  sized, 
and  the  technology  involved  is  within  the  state  of  the  art.  There  are  a num- 
ber of  speckle  theory  problems  which  need  more  detailed  examination. 
Amongst  these  we  list  the  need  for  a quantitative  understanding  of  spectral 
bandwidth  and  exposure  time  effects,  and  a theory  of  photon  shot  noise  and 
turbulence  sampling  noise  in  speckle  imagery.  The  noise  effects  are 
studied  at  a fundamental  level  in  the  next  chapter. 
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Introduction 


3.  1 

Speckle  interferometry  as  originally  demonstrated  by  1-abeyrie1 
and  speckle  imagery  as  defined  by  Knox  and  Thompson3  provide  us  with 
the  capability  of  obtaining  information  at  near  the  diffraction  limit  of  a 
telescope  when  imaging  through  turbulence.  The  key  consideration  is  that 
for  these  speckle  techniques,  while  the  spatial  frequencies  near  the  tele- 
scope's diffraction  limit  are  substantially  attenuated,  nominally  by  a factor 
of  (D/r0)a  , where  D is  the  telescope's  aperture  diameter  and  rQ  is  the 
atmospheric  turbulence  limited  coherence  diameter  3,4 , this  level  of  atten- 
uation is  not  so  great  that  it  can  not  be  compensated  in  the  speckle  proces- 
sing procedure  without  excessive  noise  effects,  if  a reasonably  large  signal- 
to-noise  ratio  is  available*.  The  key  question  thus  reduces  to  how  large  a 
signal- to -noise  ratio  is  available  in  the  speckle  process. 

Because  of  the  complex  nature  of  the  speckle  processes  the  matter 
of  evaluating  the  signal-to-noise  ratio  is  by  no  means  straightforward. 
Despite  the  existence  of  several  studies  of  noise  theory  for  speckle  inter- 
ferometry6*7 there  still  does  not  appear  to  be  an  exact  input  to  output  quanti- 
tative theory  of  noise  effects  in  the  results  of  measurements  made  using 
speckle  techniques.  In  this  chapter,  we  shall  attempt  to  develop  just  such  a 
theory  for  both  speckle  interferometry  and  speckle  imagery.  We  shall 
develop  relationships  between  the  spatial  frequency  contents  of  the  target 


* For  ordinary  long  exposure  imagery  the  spatial  frequencies  near  the  dif- 
fraction limit  are  present  in  the  recorded  image  and  at  least  nominally 
could  be  recovered  by  post  detection  compensation.  However,  in  this 
case  these  spatial  frequencies  are  attenuated  by  a factor  of  the  order  of6 
exp[-3.44  (D/r0)6/a ].  This  is,  in  general,  so  great  an  attenuation  factor 
that  no  reasonably  expected  signal-to-noise  ratio  would  be  sufficient  to 
prevent  noise  effects  from  dominating  the  post  detection  compensated  image 
if  we  attempt  to  compensate  for  attenuation  factors  as  large  as  this. 
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object,  the  number  of  detected  photons  per  short  exposure,  and  the  num- 
ber of  short  exposures,  as  our  input  quantities,  and  the  accuracy  with  which 
we  can  determine  the  magnitude  (and  phase)  of  each  spatial  frequency  com- 
ponent of  the  target  object. 


The  are  two  fundamental  sources  of  noise  (i.e.  , of  error  in  the 
final  result)  for  speckle  observations  through  turbulence.  These  two  noise 
sources  derive  from  1)  the  Poisson  statistics  of  the  photon  detection  events, 
and  from  2)  the  random  nature  of  the  turbulence  effects.  Both  of  these 
noise  sources  have  to  be  considered  in  estimating  the  accuracy  of  measure- 
ments derived  from  speckle  observations. 


In  the  next  two  sections  we  shall  develop  the  formalism  that  we  will 
use  in  this  work  to  obtain  averages  over  these  two  noise  sources.  We  shall 
then  briefly  review  speckle  interferometry  and  speckle  imagery  and  define 
the  quantities  directly  related  to  the  rms  error  in  the  final  results  of  speckle 
interferometry  and  speckle  imagery.  The  following  sections  will  be  quite 
mathematical/statistical  in  nature  and  will  develop  values  for  various 
"moments"  of  the  recorded  shot  exposure  speckle  images.  Finally,  we 
shall  pull  all  of  these  results  together  to  obtain  estimates  for  the  rms 
error  to  be  expected  in  measurement  results  extracted  from  the  speckle 
interferometry  process  and  from  the  speckle  imagery  process. 
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3.  2 Averaging  Formalism:  Turbulence 


We  shall  consider  some  target  object  which  ideally,  i.  e.  , in  the 
absence  of  turbulence  or  photon  detection  shot  noise  effects  would  produce 
an  image  pattern  ^(r)  , where  7 is  a two-dimensional  vector  denoting 
position  in  the  focal  plane  of  our  optical  system.  For  convenience  we 
shall  measure  focal  plane  position  in  terms  of  the  corresponding  field - 
of-view  position,  so  that  7 is  actually  an  angular  position  vector  measured 
in  radians. 

We  note  that  10(7)  is  everywhere  nonnegative,  i.  e.  , 

lot?)  * 0 (1) 

As  a matter  of  practical  convenience  we  define  10  (7)  to  be  so  normalized 
that 


Jd?Io0?)  = 1 (2) 

These  two  equations  allow  us,  if  we  wish  to  do  so,  to  interpret  1q  (7) 
as  a probability  density! 

Because  of  atmospheric  turbulence  effects  the  ideal  image  pattern 
will  not  be  formed.  For  the  moment,  to  limit  our  attention  to  the  random 
effects  due  to  atmospheric  turbulence  and  to  suppress  from  consideration 
the  effects  due  to  the  Poisson  statistics  of  the  photon  detection  process,  we 
shall  consider  the  target's  brightness  to  be  somehow  greatly  increased  so 
that  the  photon  detection  rate  is  virtually  infinite.  The  image  formed  in  a 
single  short  exposure  will , in  this  case,  be  random  in  nature  but  only  be- 
cause of  the  random  nature  of  the  turbulence  effects.  (With  a virtually 
infinite  photon  detection  rate  the  Poisson  statistics  of  the  photon  detection 


process  will  influence  the  image  negligibly.)  The  recorded  image  during  such 
a short  exposure,  which  we  shall  denote  by  1(7)  , will  reflect  the  random 
nature  of  the  wavefront  distortion  at  the  instant  of  that  short  exposure.  In 
exactly  the  same  way  that  we  consider  the  wavefront  distortion  to  be  a ran- 
dom function  drawn  from  a turbulence  ensemble,  we  must  similarly  consider 
I(r)  to  be  a random  function  also  drawn  from  a turbulence  ensemble. 

The  relationship  between  the  random  instantaneous  photon  noise - 
free  image  pattern,  1(7)  , and  the  ideal  image  pattern  I0(r)  is  most 
easily  expressed  in  terms  of  a random  transfer  function,  t(h)  , where 
h is  a spatial  frequency  [conjugate  to  7 , so  that  as  the  dimensions  of  7 
are  radians  (field-of-view),  the  dimensions  of  h are  radians  (cyclical) 
per  radian  lfield-of-view).  ] The  random  transfer  function  relates  the  ran- 
dom wavefront  distortion  to  the  random  image  pattern.  In  order  to  make 
use  of  this  random  transfer  function,  t(h)  , we  need  to  consider  the  spatial 
frequency  transforms  of  the  ideal  image,  I^?)  , and  of  the  random  image, 
1(7)  . We  may  write  these  as 

Mh)  = Jdr  exp(-  i7  . r)  I0(r)  , (3) 

and 

I (h)  = Jdr  exp(-  i h • r)  1(7)  , (4) 

respectively.  The  inverse  fourier  transform  relationships  are 

I0(r)  = (2n)~3  Jdn  exp  (i  h • 7)  10(>?)  , (5) 


and 


1(7)  = (2n)“3  JdJf  exp(i  k •?)  I (>?) 


(6} 
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The  random  instantaneous  transfer  function,  t(h)  , characterizing  the 
instantaneous  random  wavefront  distortion  relates  the  random  image, 

I(r)  , to  the  ideal  image  I0(r)  by  the  expression6 

I(S)  = T(K)l0(jf)  . (7) 

and  the  associated  relationship,  based  on  Eq.  's  (3)  and  (6)  in  conjunction 
with  Eq.  (7),  that 

1(f)  = (2n)”a  ^dJ?  exp(i  x •"?)  r(x)  JcT?'  exp(-i  k -r')  yr') 

= (2tt )'3  JJd^  d?'  expLi  >?•(?  - ?')]  t(T?)  I0C?')  . (8) 

For  our  purposes  Eq.  (7)  rather  than  Eq.  (8)  may  be  taken  as  the  basic 
relationship. 

In  accordance  with  standard  procedures  the  random  transfer  func- 
tion is  normalized  to  have  a value  of  unity  at  zero  spatial  frequency,  i.  e.  , 


so  that  the  random  transfer  function  can  be  written  in  terms  of  the  random 
pupil  function,  PC?)  , as9 

t(k)  = A'1  Jd?  P*  Cp  +XH)  PCp)  . <10) 

where  “p  is  a two-dimensional  vector  denoting  position  on  the  telescope  s 
entrance  aperture  (measured  in  units  of  length),  and  PCp)  represents  the 
random  pupil  function.  The  optical  wavelength  is  denoted  by  X,  so  that 
since  Tt  has  dimensions  of  radians  (cyclical)  per  radian  (field-of-view), 


Xh  is  a position  vector  with  dimensions  comparable  to  ? . The  quantity 
A in  Eq.  (9)  denotes  the  effective  optical  power  input  to  the  aperture,  i.  e. , 

A = j*d?  |P(U)|2  . (11) 

which  in  conjunction  with  Eq.  (10)  insures  the  validity  of  Eq.  (9)  . 

For  our  purposes  we  may  consider  the  optics,  per  se  , to  be 
aberration  free.  Thus,  if  we  use  the  notation  f (?)  to  denote  the  random 
complex  phase  shift  associated  with  the  turbulence  induced  random  wave  - 
front  distortion  of  some  instant,  then  we  can  write  the  pupil  function  as 

P(p)  = W(?)  exp[i  * (?)]  . (12) 

Here  W(?)  is  a function  which  serves  to  define  the  extent  of  the  tele- 
scope's entrance  aperture.  It  is  defined  by  the  equation 

1 , if  ? is  inside  the  aperture 

0 , if  ? is  outside  the  aperture  . (13) 

We  note  here  that  it  is  the  random  nature  of  ijrCp)  that  causes  the 
pupil  function,  P(?)  , to  be  a random  function.  The  random  nature  of  P(?) 
causes  the  transfer  function,  t(J?)  , to  be  a random  function.  The  random 
nature  of  t(k)  , in  its  turn,  causes  the  photon  noise-free  instantaneous 
image  pattern.  If?)  , and  its  fourier  transform,  1(h)  , to  be  random  func- 
tions. All  of  these  random  functions  are  random  in  the  sense  of  a selection 
over  the  ensemble  of  turbulence  realizations.  We  shall  use  the  notation 
(f)T  to  denote  an  average  of  some  random  function,  f , over  turbulence, 
i.  e.  , over  the  ensemble  of  all  possible  realizations  of  the  turbulence  pat- 
tern and  the  associated  values  of  f . 


W<P>  = 
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We  would  write,  for  example,  for  the  mean  and  mean- square 
optical  transfer  function  at  spatial  frequency  >?,  (t(k))t  , and  (|t(>?)|  )T  , 

respectively.  Making  use  of  Eq.  (10),  and  of  the  fact  that  the  ensemble 
averaging  process  is  linear  and  therefore  cominutes  with  the  integration 
process,  we  could  write 

<t(h))t  = <A”1  J*d?  P*  C?  + X*)  P(?)>T 

= A-i  Jd?  <P*<?  + \k)  P(?)>t  . (14) 

Similarly  we  could  write 

< | t(h)|  8 >T  = <CA-»  J*d?  P*(?  + X*)  P(?)]* 

X LA-i  P*(p'  + X*)  PC?/)]>t 

= <A*«  J*J*d?  d?'  PCp  + X*)  P*(p) 

X P*  tf'+X*)  P(?)>T 

= A-aJJd^d^  <P(*p  + X#)  P*(p)  P*(p/+X^)P(7f,)>T  . (15) 

The  farther  reduction  of  expressions  like  those  in  the  right-hand- side  of 
Eq.  's  (14)  and  (15)  is  a matter  which  we  shall  consider  later  in  discussing 
the  averaging  over  turbulence.  At  this  point,  we  may  consider  our  formalism 
for  separating  the  random  turbulence  induced  effects  from  the  random  photon 
detection  effects  as  established,  as  well  as  our  formalism  for  writing  the 
turbulence  average.  We  are  now  ready  to  turn  to  the  development  of  a pro- 
cedure for  handling  the  random  nature  of  the  photon  detection  process.  We 
take  this  up  in  the  next  section. 
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3.  3 Averaging  Formalism:  Photon  Detection 


We  note  here  that  the  random  photon  noise-free  instantaneous  image 
pattern.  I("r)  , since  it  corresponds  to  a local  image  density  is  everywhere 
nonnegative,  i.  e.  , 

U?)  * 0 . (16) 

Moreover,  if  we  make  use  of  Eq.  's  (3),  (6),  and  17),  we  can  write 

J*d?  1(7)  = (2n)“a  J'J'Jdn  d"r  d?#  exp[-  in  *(r#-7)] 

X T (h)  10(?')  . (17) 

This  expression  can  be  reduced  by  making  use  of  the  fact  that  the  repeated 
fourier  integral  recovers  the  starting  function.  With  5 , , and  y'  as  n- 

dimensional  vectors  we  can  write  this  property  of  the  repeated  fourier 
integral  as 


J’J'dS  dy  expfcfc  i x • (?  - y')]  h(y')  = (2tt)»  h(y)  , (18) 

where  h(y)  is  any  reasonably  well  behaved  function.  Choosing  the  minus 
part  of  the  plus/minus- sign  in  Eq.  (18),  associating  h in  Eq.  (17)  with 
in  Eq.  (18),  ? in  Eq.  (17)  with  x in  Eq.  (18)  , letting  ^ in  Eq.  (18)  be 
set  equal  to  zero,  and  equating  h(7')  in  Eq.  (18)  with  t(h)  I0(r/)  in  Eq. 
(17),  we  can  use  Eq.  (18)  to  allow  us  to,  in  effect,  perform  the  h-  and 
7 -integrations  in  Eq.  (17).  Thus  we  can  rewrite  Eq.  (17)  as 

/dfltf)  = (2TT)-»JV?'t(2TT)3  T(0)lo(?')} 

= T(0)Jd*?'  10(7') 
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(19) 


Making  use  of  Eq. ' s (2)  and  (9)  this  result  can  be  reduced  to  the  form 

J*dr  I(T)  = 1 . (20) 

This  result  is  very  significant  since  it,  in  conjunction  with  Eq.  (16),  allows 
us  to  interpret  the  photon  noise-free  instantaneous  image  pattern,  I (?)  , 
as  a probability  density!  1 (?)  is,  in  fact,  the  probability  density  govern- 
ing the  position  distribution  for  each  detected  photon.  We  shall  use  it  to 
define  the  probability  of  a photon  being  detected  at  various  positions  in  the 
focal  plane. 

The  random  aspects  of  the  photon  detection  process  can  be  usefully 
separated  into  two  parts.  The  first  part  concerns  the  number  of  photons 
detected  in  forming  a single  short  exposure.  Through  we  expect  this  num- 
ber to  have  the  mean  value  N , the  actual  quantity  detected  in  any  one 
short  exposure,  N , is  a random  number.  The  second  part  concerns  the 
location  on  the  focal  plane  where  each  of  the  N photon  detection  events 
actually  occurs.  At  least  nominally  we  expect  these  events  to  be  distri- 
buted in  a way  that  matches  the  image  pattern/probability  density,  I(?)  . 

In  fact,  the  distribution  on  the  focal  plane  will  be  random  and  the  match 
to  the  image  pattern,  1(?)  , will  be  imperfect.  In  treating  photon  detection 
we  shall  have  to  take  account  of  these  two  things  jointly  so  as  to  fully  account 
for  the  random  nature  of  the  photon  detection  process. 

In  treating  the  random  nature  of  the  photon  detection  process  we 
shall,  in  general,  be  treating  some  random  function  ^(?)  which  defines 
the  distribution  of  photon  detection  events  on  the  focal  plane  for  a single 
short  exposure.  We  can  write 

N 

*<?)  = y «(?  - *?,)  . (2D 

L- 

t*  1 
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where  6(r  - ft)  is  the  Dirac  delta  function.  Here  N is  a random  num- 
ber representing  the  number  of  photon  detection  events  in  the  short  expo- 
sure, and  rt  is  a random  position  vector  denoting  the  location  of  the  i™ 
photon  detection  event.  The  value  of  N is  governed  by  a Poisson  distri- 
bution with  mean  value  TC  , i.  e.  , 

Prob(N)  = —j  exp(-TT)  , (22) 

while  the  distribution  of  rt  is  governed  by  the  image  pattern/probability 
density  I (r,)  . It  follows  that  the  mean  value  of  some  function  of  N , 
e.g.  , f(N)  , averaged  over  the  photon  detection  statistics  is 

<f(N)>,  = 7 “ exp(-N)  f(N)  , (23) 

NaO 


while  the  average  of  some  function  of  rt  , e.g. , g(7,)  , averaged  over 
the  photon  detection  statistics  is 

<gC?t)>,  = I(tt)  g(?t) 

= Jdr  I(r)  g(T)  . (24) 

In  both  Eq.  (23)  and  Eq.  (24)  the  angle  bracket  notation  subscripted  P , 
i.  e.  , ( )p  , denotes  an  average  over  the  photon  statistics. 

It  is  particularly  significant  to  note  that  the  photon  detection  pro- 
cess is  such  that  each  detection  event  is  independent  of  each  other  detection 
event.  Thus,  with  i ^ i'  , for  some  function  G (?t  , r^)  we  can  write 

«G(?t  .■?„)>  = /Id*?,  d?t,  Iff,)  I(*„)  G(?t  , ?,,)  , (25) 
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and  if  C(?t,  r„)  is  factorable,  i.  e.  , we  can  write 


Gtfj.'t#)  = gir?t)  g8(?t/)  . (26) 

then 

<c(?t . ?u))  = ;;d?4  dsl#  lc?,)  K?t/)  8a(?1#) 

= C ;<*»  Krt)  gl(^t)}  c,fd?l#  k?w)  ga(?1#)3 

= { J*d?  I(?)  g ,(■?))  { Jd?  I(r)  g ,(■?)) 

= <gx(r)>„  <gaC?)>p  . (27) 

In  our  analysis  we  shall  be  interested  in  the  average  over  photon 
detection  statistics  of  the  produce  of  two  functions  of  the  same  randomly 
chosen  short  exposure.  The  functions  will  be  of  the  form 

N 

Gi  = £ 8i(rt)  , (28a) 

1*1 

N 

Ga  = Y g8C?,)  . (28b) 

i»i 

where  N has  the  same  random  value  in  both  expressions,  and  the  random 
vectors  , ?a  , ?3  , . . . , rN  are  identical  in  the  two  expressions.  We 
shall  be  interested  in  evaluating  the  quantity 

N N 

Ga)p  = ^ ^ 81^1)  y 8a(^i/))p 

1 = 1 i/»i 
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It  is  convenient  to  separate  the  N®-  terms  in  this  double  sum  into 
N(N-l)  - terms  for  which  i ^ i'  and  the  N-terms  for  which  i = i'  . 
Thus  we  can  write 

N 

<Gi  Ga>P  = < y gi^t)  ga('?1/)>f 


+ < J 8a  • 


Making  use  of  Eq.  (24)  we  can  see  that  all  N-terms  in  the  second  en- 
semble average  in  Eq.  (30)  reduce  to  the  same  form  so  that  we  can  write 

<GxGa  >p  = < T 8xC?t)  ga f?w)>, 

/ — i 

1,1/al 

i^i/ 

+ <N>p  ;d?I(?)gl(?)ga(T»)  . (31) 

Making  use  of  Eq.  (27)  we  can  see  that  when  i / i*  , the  photon  detection 
event  average  over  gjf?^  ga(rt/)  in  Eq.  (30)  can  be  separated  into  two 
distinct  averages.  Thus,  the  N(N-l)-terms  in  this  double  sum  can  be 
written  in  such  a way  that  we  get  in  place  of  Eq.  (31), 

<Gx  Ga )p  = <N(N-1)>„  [ J*d?  1(f)  gxC?)3£ /d?  1(f)  ga (?) } 


+ <N>p£j*d?lC?)  gl(?)  ga(?))  . 
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In  developing  Eq.  's  (31)  and  (32)  we  have  taken  advantage  of  the  fact  that 
the  integral  over  r , i.  e. , Jd"?  I(?)  ....  is,  in  effect,  an  average  over 
the  photon  detection  event  statistics  in  and  of  itself,  and  therefore,  since 
it  is  not  a random  quantity,  need  not  be  retained  inside  the  angle  brackets. 
Only  the  quantities  N and  N(N-l)  , which  count  the  number  of  terms  in- 
volved in  the  sum,  remain  inside  the  angle  brackets.  These  quantities  can 
be  evaluated  making  use  of  Eq.  (23)  . In  fact,  we  find  using  Eq.  (23)  that 

<N>„  = N . (33a) 

<N(N-  l))p  = N3  , (33b) 

and  that  in  general,  the  QTh  factorial  moment  of  N,  i.  e.  , (N(N-l)  (N-2) 

. . . (N-Q+l))p  has  a value  of 

<N(N-1)  (N-2)  . . . (N-Q+l)>„  = N9  . (34) 

Making  use  of  these  techniques  suitably  extended  to  allow  for  the  evaluation 
of  the  mean  value  of  the  produce  of  not  just  two  functions  of  the  form  of  Gx 
and  Ga  , but  of  the  produce  of  up  to  four  such  functions,  we  shall  be  able 
to  develop  the  photon  detection  statistics  for  all  of  the  quantities  of  interest 
to  us  in  speckle  interferometry  and  speckle  imagery.  We  carry  out  the 
basic  evaluations  in  Section  5.  Before  turning  to  that,  however,  we  shall 
first  briefly  comment  on  the  formalism  for  the  joint  photon-detection-event 
and  wavefront-distortion  statistics.  We  take  this  up  in  the  next  section. 
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Averaging  Formaliam:  Photon  Detection  and  Wavefront  Distortion 


In  practice  we  shall  be  interested  in  averages  taken  over  the  statis~ 
tical  nature  of  both  the  wavefront  distortion  process  and  the  photon  detection 
process.  We  shall  use  the  notation  of  unsubscripted  angle  brackets,  vis.  , 

( ) , to  denote  such  a joint  average.  For  evaluation  we  shall  understand 
that  we  may  replace  this  average  notation  with  the  notation  <(  ) \ , with  the 
understanding  that  the  photon  statistics  average  is  always  taken  within  the 
turbulence  (i.  e.  , wavefront  distortion)  average.  This  is  because  the  photon 
statistics  average  may  be  expected  to  be  expressed  in  terms  of  I(?)  , the 
image  pattern/ probability  density.  However,  this  image  pattern  is  random 
in  the  sense  of  its  dependence  on  turbulence  and  so  the  photon  detection 
average  taken  over  I(r)  , must  then  be  further  averaged  to  account  for  the 
random  variability  of  I(r)  associated  with  turbulence  effects.  Hence  the 
requirement  that  the  photon  detection  average,  ( )p  , be  enclosed  within 
the  angle  brackets,  ( )T  , denoting  the  average  over  turbulence  effects. 

With  this  matter  settled,  we  are  now  ready  to  turn  our  attention  to 
evaluating  some  of  the  photon  detection  statistical  averages  that  we  will 
need  in  evaluating  the  speckle  statistics.  We  take  this  up  in  the  next 
section. 


Spatial  Frequency  Spectral  Moments:  Photon  Detection  Averaging 


3.  5 

Speckle  interferometry  and  speckle  imagery  are  both  based  on  ma- 
nipulation of  the  spatial  frequency  transform  of  the  short  exposure  image. 

It  will  therefore  prove  convenient  in  our  subsequent  work  if  we  develop 
here  a set  of  expressions  for  the  photon  detection  event  statistical  averages 
of  various  moments  of  the  random  amplitudes  associated  with  the  fourier 
transform  of  a random  short  exposure  image.  We  develop  these  results  in 
the  following. 

The  basic  random  short  exposure  image  pattern  is  represented  by 
A?(?)  , which  was  defined  by  Eq.  (21).  In  that  equation  N , as  we  recall, 
denotes  the  random  number  of  photon  detection  events  making  up  the  image, 
and  r,  denotes  the  location  of  the  i™  photon  detection  event.  The  spatial 
distribution  of  ?,  is  governed  by  the  probability  density  H^)  , which  so 
long  as  we  are  considering  only  the  photon  detection  statistics  may  be 
considered  to  be  a well  defined  quantity.  [It  is  only  when  we  turn  to  con- 
sideration of  averages  over  turbulence  that  we  have  to  take  account  of  the 
fact  that  I(?)  is  a random  function.  In  the  work  of  this  section  we  may 
ignore  that  fact.  ] The  fourier  transform  of  the  random  pattern  at  spatial 
frequency  may  be  written  as, 

R(^j)  = Jdr  expl-iKj-r)  ^(r)  . (35) 

When  we  substitute  Eq.  (21)  into  Eq.  (35)  , interchange  the  order  of  inte- 
gration and  summation,  and  then  use  the  well  known  properties  of  the 
Dirac  delta  function  to  allow  us  to  carry  out  the  integration,  we  obtain 
the  result  that 

N 

R(J?j)  = ^ exp(-iKj«r,)  . (3b) 

t-  l 
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We  shall  be  interested  in  the  average  over  photon  detection  statistics  of 
R(h4)  and  of  the  produce  of  such  quantities  for  different  values  of  . 


Because  of  the  fact  that 

. (37) 

we  can  restrict  our  attention  to  the  following  fo.i  averages;  namely 

= <R(*i)>,  . (38a) 

Ma(*x.Ha)  = <R(5fx)  R(Ha)  R(Sa)>,  . (38b) 

Hs)  = (r(hj)  R(*a)  R(^a))^  » (38c) 

M*(*i  , J?a  , k3  , h4)  = <R (ifj)  R(h8)  R(*s)  R(k4)>,  . (38d) 


Storting  with  the  evaluation  of  Mx  , substituting  Eq.  (36)  into  Eq.  (38a) 
we  get 

N 

Mi  fl*»)  = < 7 exp(-  i?x  . ?t)>,  . (39) 

Making  use  of  Eq.  (24)  and  the  formalism  developed  in  Section  3,  we  can 
rewrite  this  as 

Mx(Tfx)  = <N>,  J* d?  exp(-  i hx  • t)  1(?) 

= NI^j)  , (40) 


where 


1(h)  = Jd?  exp(-  i h •?)  If?)  . (41) 

We  note  that  just  as  If?)  is  a randomly  selected  function,  held  constant 
within  the  photon  detection  statistics,  similarly  I(?)  is  also  randomly 
selected  function  held  constant  within  the  photon  detection  statistics. 

[Both  I(r)  and  1(h)  are  randomly  selected  from  the  turbulence  ensemble.] 
Just  as  If?)  corresponds  to  a probability  density,  1(h)  represents  the 
associated  generating  function. 

The  evaluation  of  Ma  follows  almostly  directly  from  Eq.  (32) 
and  the  immediately  preceding  type  of  manipulations.  Starting  with  Eq. 
(38b)  and  substituting  Eq.  (36)  twice,  we  can  write 

N N 

Ma(>?i,Ha)  = < \ exp(-  i Jfj  • rt)  V exp(-  i h8  • r,,)), 

i=x  i'=i 

N 

= < ^ exp(-  i Hx  *?t)  exp(-  i Ha  •■?!#)>,  . (42) 

t,i'=  l 

Following  the  procedure  demonstrated  in  developing  Eq.  (32)  from  Eq. 

(29)  , we  start  by  separating  the  double  sum  in  Eq.  (42)  into  the  N(N-l)  - 
terms  for  which  i ^ i#  and  the  N-terms  for  which  i = i#  . Thus  we  can 
write 

N 

M»(*i**a)  = < Y exp(-  i ^‘Tj)  exp(-  i Ha  •?,/)>, 

i,i/=  l 
i#i/ 


46 


: ~ | ~ ^ i a imi  i 


(43) 


Proceeding  just  as  in  the  development  of  Eq.  (32)  we  can  replace  each  of 
these  exponentials  by  their  average  over  the  probability  density  If?j)  or 
I(?j/)  • and  take  the  average  value  outside  of  the  angle  brackets.  The  two 
sums,  each  with  unity  as  the  summand,  reduce  to  N(N-l)  and  N . Thus, 
Eq.  (43)  can  be  reduced  to  the  form 

Ma(Hx  . *a)  = <N(N-1)X  [ J*drt  expt-i^-?,)  I(?t)3 

X { Jd?5/exp(-i  h3  •?„)  I (?,,)] 

+ <N)P  C Jd?,  exp[-i  (^  + h8)  • rj  I(rt)3  . (44) 

Making  use  of  Eq.  's  (34)  and  (41)  we  can  reduce  this  to 

MjtHj.Hj)  = I(5TX)  f(H8)  +K  ItHj+Ha)  . (45) 

This  same  procedure  can  now  be  applied  for  the  evaluation  of  the  higher 
moments,  Ma  and  M4  . 

In  evaluating  Ma  we  will  encounter  a triple  sum  on  i,  i*,  and  in 
and  will  have  to  consider  a more  complex  pattern  of  subdivision  of  the  sum. 
If  we  substitute  Eq.  (36)  into  Eq.  (38c)  we  get 

N N 

MjfKj  , Ha  , h3)  = < y exp(-i  -tj)  v exp(-i  ks  • rt/) 

i*l  i/=  l 

N 

x Y exp(-iH3  , (46) 

i*=  l 
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which  when  we  form  a triple  sum  from  the  product  of  three  sums,  takes 
the  form 


! 


N 

Mg (hj  i Kg  i Jfj)  ~ ^ J exp  ( - i Hj  • ^j)  exP( ” t 

x exp  (-i  k3  'ri0))f  , (47) 

In  this  case  we  separate  the  triple  sum,  consisting  of  N3 -terms,  into  a 
sum  on  the  N(N-  1)  (N-2)-terms  [equal  to  (N3  -3N3  + 2N)-terms]  for  which 
[i  ji  i1  iw]  , three  separate  double  sums  of  N(N-l)-terms  [equal  to  a 
total  of  (3N2 -3N)-terms]  for  which  [i  = i-  ] , or  [i  = i0^  if  ] , or  (i#  = i0 4 i]. 
and  the  single  sum  of  N-terms  for  which  {i  = i'  = i0  ].  Thus,  we  can  write 

N 

Mg^  , h8  , s?3)  = < ) exp(-i  Hx*  Tj)  exp(-i  Ka- rt/)  exp(-i  Kg-?l#f))p 

1,1', 1**  1 

1 * I'M  l" 


+ < y expt-i^  + Ka)  • r,]  exp(-i  h3  •■?„)>, 


i,  i*=  x 

irf  i* 


+ < Y exp[-i(H1  + K3)-r1]  exp(-i  h8- tj,)), 

i,i'«  x 

1 4 1' 

N 

+ < y exp[-i(Kg+ Kg)  • rt/]  exp(-i  Hx  • rt))p 


i,i'*  l 

l*  i' 
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(48) 


I* 

•«( 

n 

f 

i. 


Averaging  over  the  portion  of  the  photon  detection  statistics  which  concerns 
the  location  at  which  each  of  the  N-photons  is  detected,  which  is  governed 
by  the  probability  density  l(r)  , and  using  the  fact  that  the  random  variables 
r,  , 7„  , and  7,#  as  they  appear  in  each  angle  bracket  term  in  Eq.  (48) 
are  distributed  independently  of  the  others,  we  can  write  in  place  of  Eq.  (48) 
that 


Ma(Tf1.  h8,  *3)  = <N(N-l)(N-2)>,  {/d7,  exp(-i  ?,.?,)  1(7,)) 

* [J’<1?1/exp(-i»?8.?1#)  1(7,,))  {/dr„  exp(-i  h, -7,„)  1 (?,,)} 

+ <N(N-1)V  Cjd^expMOf,  + *„).*,]  1(7,)) 

X «xp(-i  S3  -7,#)  1(7,,)} 

+ <N(N-l)>p  {/dr,  exp[-i(?x  + *a)-7,]  1(7,)) 

X </<*?* / e*P(-i  h8-?„)  1(7,,)) 

+ (N(N-l))p  {/dr,,  expC-i  (78 + ?,)•?„ ] 1(7,,)) 

X {/d*?,  oxp(-i  Hx  •?,)  1(7,)) 

+ <N>,  t/dt,  exp[-i(^  + ^8  + ^s).?,]  1(7,))  . (49) 

Now,  if  we  make  use  of  Eq.  's  (34)  and  (41)  we  can  reduce  this  to  the  much 
more  compact  result  that 

M.Of,  , J?8  , 7ta)  = N3  l(^)  I(7t)  UZa) 
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(50) 


+ N 8 1(HX  + Ha)  I (Ha)  + N 3 I(hx  + 73)  I(Sa) 


+ T7s  I(h3  + K3)  I (»?x)  + N I(HX+Ha+Ha) 

This  same  sort  of  procedure  now  has  to  be  carried  out  to  evaluate  the 
fourth  moment,  M4  . 

In  evaluating  M4  we  initially  form  a product  of  four  sums  which 
we  rewrite  as  a quadruple  sum  of  N4  - terms.  Thus,  we  have 

M4(J?x,  >?a,  >t3.  h4)  = < exp(-iHl-r,)  V exp( -i  Ha  • r„) 

1 m 1 i'=  1 

N 

X ^ exp(-i  Jt3 -1lM)  y exp(-i  h4 

L-  i— 

Hsl  1*=1 


= < exp(-i  •?,)  exp(-i  Ha.f1#) 

1,  W,l#,l»=  1 

X exp(-i  J?3  . (51) 

This  sum  of  N 4 -terms  can  be  separated  into  the  following  subsidiary  sums: 

1)  a sum  of  N(N-l)  (N-2)  (N-3)-terms  for  which  {i  ji  i*  jl  iV  i*},  2)  six  sums 
of  N(N-l)  (N-2)-terms  for  which  (i  = i'^  i,#^i*f),  or  [i  = i*^  iV  i*},  or 
{i  = i*7  iVi')  , or  [i'  = i i ? i").  or  [i'  = i*7  i / i"),  or  {i#=  i *7  i 4 i'}, 

3)  three  sums  of  N(N-  l)-terms  for  which  { i = i#  ^ i*  = i*  ) , or  { i = i*  ^ i'  = i*  } , 

or  [i  = i*7  i*  = i*},  4)  four  sums  of  N (N  - 1)  - terms  for 

which  (i  = i'  = i*  f i*  ],  or  {i  = i'  = i*  ^ \*  ) , or  [i  = i*  = im  ^ i*  ) , or 

Ji*  = i*  = \m  £ i},  and  5)  one  sum  of  N- terms  for  which  {i  = i*  = i#  = i**  ] . 

Following  exactly  the  same  procedures  as  we  used  in  going  from  Eq.  (47) 

- 50  - 


to  Eq.  (50)  this  separation  of  the  terms  in  the  sum  in  Eq.  (51)  and  the 
subsequent  reduction  will  yield  the  result  that 


m4  (kx  . k#  . h3  . **)  = n4  M^)  r (>g  r(Tf3)  r(»4) 

+ tt3  i(5?x  + >{a)  i(>?3)  i(j?4)  + T73  + jt,)  i(j?8)  ro?s) 

+ N3  l(Xx  + h4)  I(JTa)  I(m*3)  + N 3 I (h8  + *s)  f(^)  f(*4) 


+ N3  f(H8  + J?4)  I(J?X)  I(*3)  + N3  I(*3  + J?4)  I(KX)  l(J?a) 

+ N®  1(»?1  + J?a)  1(J?3  + JT4)  + N®  Itfj  + h3)  I(k8  + J?4) 

+ K3  + >l4)  1 (h8  + Jig)  + N3  I (H1  + Jig  + Jig)  1 (Jl4) 

+ TT3  HKj  + Jlp  + K4)  l" (Jig)  + TT3  1 (Hj  + Jig  + Jl4)  l(Jl8) 

+ Na  f(j?a  + Kg  + J?4)  I (Jfj  + N I (Jfx  + JTa  + Hg  + H4)  . (52) 

With  these  results  for  the  first  four  spectral  moments  in  hand,  i.  e.  , M1  , 
Ma,  Mg  , and  M4  , as  given  by  Eq.  's  (40),  (45),  (50),  and  (52)  we  are  now 
ready  to  consider  the  quantities  calculated  in  speckle  interferometry  and 
speckle  imagery,  as  averaged  over  the  photon  detection  statistics.  This 
is  treated  in  the  next  section. 
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3.6 


Unbiased  Speckle  Measurement  Quantities 


In  speckle  interferometry  and  speckle  imagery  the  basic  mea- 
surement is  the  fourier  transform  at  spatial  frequency  h , of  the  random 
image  pattern  recorded  during  a single  short  exposure.  This  corresponds 
to  the  quantity  R(>t)  defined  in  Eq.  (35).  Data  reduction  is  based  on  cal- 
culating the  mean  value  of  a quantity  like  R*(h)  R(h)  for  speckle  interfero- 
metry and  of  a quantity  like  R*(h7)R(x)  for  speckle  imagery.  The  averag- 
ing is  over  a series  of  short  exposures  and  thus  is,  in  effect,  over  both 
photon  detection  statistics  and  over  turbulence  statistics.  Since  the  expres- 
sion R*(>?7)  R(>?)  for  speckle  imagery  obviously  subsumes  the  expression 
R*(x)  R(x)  for  speckle  interferometry  we  shall  treat  only  the  former  quantity. 

The  average  over  many  short  exposures  for  this  quantity,  i.  e.  , 

(R*(>t7)  R(>?))  , which  represents  an  average  over  both  the  turbulence  de- 
pendence and  the  photon  detection  effects,  can  be  rewritten  using  the  for- 
malism of  Section  4,  as 

(RV)R|k))  = < <R*(hO  R(h)),)t  . (53) 

Making  use  of  Eq.  's  (37),  (38b),  and  (45)  we  can  rewrite  this  as 

(R*(5?7)  R(k)>  = <T431(-h')  1(h)  +T7f(J?  - h7)>t  . (54) 

Now  making  use  of  Eq.  (7)  and  removing  from  inside  the  angle  brackets 
all  the  terms  which  do  not  depend  in  a statistical  way  on  the  wavefront 
distortion,  we  can  rewrite  Eq.  (54)  as 

<R*(S')  R(h))  = N8  Ie(-S')  I0(S)  <t( -k7)  t(k))t 

+ NI0(J?  -h)  <T(J?  -*')>T  . (55) 
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We  note  from  Eq.  's  (4)  and  (10)  that 


i0(-h')  = r0*(*) 

and  that 


(56) 


t(-k')  = , (57) 

so  that  Eq.  (55)  can  be  rewritten  as 

<R*on  R(*)>  = Naf0 *(;?')  I0(S)  <T*£')  t(S))t 

+ N10(h-5T)  . (58) 

This  result  is,  in  a sense,  potentially  embarrassing  in  as  much  as  it 
indicates  that  R*^)  R(h)  is  a "biased"  quantity.  We  note  that  10  (h) 
represents  the  spatial  frequency  spectral  amplitude,  at  spatial  frequency 
h , of  the  object  whose  pattern,  I0(?)  , we  are  trying  to  determine.  We 
are  entitled  to  expect  that  if  the  object  has  no  spectral  amplitude  at  spatial 
frequency  k or  >?'  then  the  average  we  are  evaluating,  nominally 
(R*(>?#)  R(h)>  , should  have  zero  value.  However,  we  can  see  from  con- 
sideration of  Eq.  (58)  that  because  of  shot  noise  effects,  this  average  does 
not  vanish  when  ^(h')  or  1^(5?)  is  zero.  There  is  a residual  bias  term 
of  the  order  of  N Ip  (5?  - k')  (t(h-h/))t  . Even  with  large  values  of  N , 
because  of  the  disparity  between  <t(h->?/))t  and  (t*  (l?7)  t(J?))t  when 
h and  h*  are  nearly  equal  and  h is  large  (i.  e.  , near  the  telescope  dif- 
fraction limit  cutoff), this  bias  term  can  be  substantial  compared  to  the 
"principle"  term  in  Eq.  (58). 
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In  order  to  avoid  this  bias  term  the  speckle  processor  should 
evaluate  the  average  over  a series  of  short  exposures  of 

S (J?  , >?')  = R*0?')  R(S)  - R(*  -Jf')  . (59) 

It  is  easy  to  show  that  the  average  of  this  quantity  leads  to  an  unbiased 
result  of  the  basic  form  we  expect.  We  can  write,  using  the  formalism 
of  Section  4,  that 

<S(h.  *')>  = < <S(Jf , H#))p  )T  . (60) 

Substituting  E{.  (59)  into  Eq.  (60)  we  get 

<S(h.  *')>  = <<R*(x')  R(x)),X  - <<R(>f  - S'))p>T  . (61) 

Making  use  of  Eq.  (58)  as  well  as  of  Eq.  's  (38a)  and  (40)  allows  us  to 

rewrite  Eq.  (61)  as 

<S(S.  >f')>  = Na  I0*(S')  I0(*)  <T*<*')  t(h))t 

+ Kr0 (*-£')  <t(S  - >T')>t 

“ ^ <f(x  - h/)>t  . (62) 

Now  making  use  of  Eq.  (7)  to  allow  us  to  write 

<r<j?  - ^)>T  = r0o?-*')  <t(j?  - ^)>T  , (63) 

we  can  reduce  Eq.  (62)  to  the  form 

<S(if,  S')>  = i!s  ro(3)  <T*(lf')  t(k)>,  . (64) 


[ j 

jj 

1 
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This  clearly  is  an  unbiased  average  containing  exactly  the  form  of  de- 
pendence we  expected,  i.  e.  , a dependence  on  the  objects  spectral  ampli- 
tudes at  spatial  frequencies  >?  and  Jf#  , and  on  the  turbulence  averaged 
product  of  the  transfer  functions  at  these  two  spatial  frequencies. 

We  are  now  ready  to  ask  how  accurately  does  the  average  of  S 
over  a finite  series  of  short  exposures  approximate  the  expected  value. 
We  expect  deviations  due  to  the  fact  that  we  have  only  a finite  number  of 
short  exposures,  thus  providing  an  imperfect  average  over  turbulence, 
and  due  to  the  fact  that  there  are  only  a finite  number  of  short  exposures 
and  a finite  number  of  photon  detection  events  per  short  exposure,  thus 
providing  an  imperfect  average  over  the  statistics  of  the  photon  detection 
process.  These  matters  are  taken  up  in  the  next  two  sections. 
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3.  7 Speckle  Interferometry  Measurement  Accuracy 

For  evaluation  of  measurement  accuracy  in  speckle  inteferometry 
we  shall  be  concerned  with  the  quantity 

CJh.Z')  = <[S(>?\S')  - <S<*',S')>]*  [S(S?,*)  - (S(>?  , k)>]) 

= <S*(J?',>?')  S(>? , #))  - *?'))  <S(h  , h)>  . (65) 

In  Section  9 we  shall  show  exactly  how  this  quantity,  CL  (h  , h')  , is 
related  to  the  accuracy  of  a speckle  interferometry  and  will  carry  out 
a sample  evaluation.  In  this  section  we  shall  be  concerned  with  evaluation 
°f  CL  (h  , 7t*)  averaged  over  photon  shot  noise. 

We  start  our  analysis  with  consideration  of  the  quantity 
* ** ) S(h  , k))  • Making  use  of  Eq.  (59)  we  can  write  this  as 

<S*(^,^)  S(h  .*)>  = <[R*(S')  Rtf')  - R*(0)] 

X [R*  (*)  Rtf)  - R(0)]>  . (66) 

It  follows  from  consideration  of  Eq.  's  (21)  and  (35)  [or  from  considera- 
tion of  Eq.  (37)],  that 

R*(0)  = R(0)  . (67) 

This  allows  us  to  rewrite  Eq.  (66)  in  the  form 

<S*(JT,^)  SO?  , *)>  = <R*(h')  R (5?')  R*(J?)  R (>?)) 

- <R*(*')  Rtf')  R(0)>  - <R*tf)  Rtf)  R(0)> 
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+ <R(0)  H(0)>  . 


(68) 


^ ^ccordincc  with  the  not^tioncil  provisions  developed  in  Section  4,  we 
may  make  explicit  the  distinction  between  the  turbulence  average  and  the 
photon  average  in  Eq.  (68)  . Thus,  we  write 

<S*(5?'  , *')  S(Z , *)>  = «R*(7f')  R(>?')  R*(J? ) R(5?))p  >T 

- «R*(*')  R(JT)  R(0)>p>T  - <<R*(>?)  R(5?)  R(0))p)t 

+ (<R(0)  R(0)>p)T  . (69) 

Making  use  of  Eq.  !s  (37)  and  (38b,  c,  and  d)  we  can  rewrite  this  as 
<S*(J?',S/)  S(*,;?)>  = nm*(-S',*',-  TT,  *)>T 

- <M3(-k'  , 0))T  - <M3(.Jf,  it,  0)X 


+ <Mg(0, 0)>  . (70) 

Substituting  Eq.  's  (45),  (50),  and  (52)  into  Eq.  (70)  we  get 
<S *(*',*')  S (it.it))  = N4<1(-h0  I(hO  1(-J?)  I(*)>t 

+ T73{<1(0)  l(-it)  i(>?)>T  + <r (-*'-*)  W)  I(h)>t 
+ <I(->?'+*)  1(5?')  T(-7)>t  + <X(5?'-5?)  I( -5?')  f(^)>T 
+ <I(ff'  + it)  X(-V')  I(-h)>t  + <I(0)  f(-^)  1(^)>T  ] 

+ N2  { <f(0)  f(0))T  + <!'(->?'-;?)  f(jr+JT)>T 
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+ (n-*?*  + x)  f(jr-  h)>t  + <i(-^)  r(j?)>T 

+ <I(*)  f(^)>T  + <f(-HO  ion>T  + <?(*')  1(-S')>T  } 

+ H <I(0)>T 

- TT3<I(.^)  1(70  I(0)>T 

- TTa  £<I<0)  I(0)>T  + <I( -7T)  I(r»T  + <I(>?')  I(-^)>T} 

- N <I(0)>T 

- H3<I(-h)  1(h)  I(0))t 

- Na{<i(o)  f(o)>T  + i(j?)>T  + <i(^)  r(-^)x } 

- 7?  <I(0)>T 

I 

+ TTa<i(0)  I(0)>T +N  <I(0))T  . (71) 

A large  number  of  the  terms  in  Eq.  (70)  cancel  so  that  the  expression  can 
be  reduced  to  the  form 


I 


<s*(tf.  H')  S(*  .*)>  = T?4  <f(-^)  W)  i(-j?)  f(S)>T 
+ TT3£<f(-H#-J?)  f(^)  i(^)>T+<f(-^+j?)  f(Jf')  i(j?))T) 

+ <r(*'->?)  f(-^)  i(S)>T  + <f(^+ >?)  I (-S')  f(-S)>T } 

+ 7?3£<I(-Z'-S)  I(^+JT)>t  + <f(-S'+S)f(S'-S))T)  . (72) 


- 58  - 


Making  use  of  Eq.  (7)  to  allow  further  reduction  of  Eq.  (72),  and 
then  combining  that  result  with  that  of  Eq.  (64),  we  can  now  rewrite  Eq. 
(65)  as 


cL(*,S')  = N4|r0(*')|a  |r0(*)n<|T(*')|a  |t(*)|3>t 

- <|T<*')|3>T<|T(*)|a>T) 

+ TT^iVHh'  +*)  I0(*')  I0(h)  <t*(^  + k)  t(S')  t(>?)>t 
+ IW-*)  I0(S')  I0*(>?)  <r*(>?'-ff)  t^)  t*(*)>t 

+ l0  (S'+k)  l0*(S')  I0*(S)  <T0?'  +1?)  T*(*')  T*(T?)>t ] 

+ N3  { 1 10  (k*  + Jt) | 3 <|t(k/+  h)|3)t 
+ |r0c?f-»)|*<|T(if'-*)|a>T}  . (73) 

With  this  result  in  hand  we  are  now  ready  to  consider  the  averages  over 
the  turbulence  statistics.  This  is  taken  up  in  the  next  section. 


3.  8 Averaging  Over  Turbulence 


For  our  purposes  here  it  will  be  sufficient  to  use  approximate 
results  for  the  turbulence  averages.  Asa  practical  matter  we  may 
further  restrict  our  attention  to  values  of  h and  n#  that  are  much 
larger  than  rQ/X.  . Accordingly,  we  may  make  use  of  the  approximation^10 

<|t(h)|3>t  « 0.435  (r0/D)3TDL(5?)  , (74) 

where  r0  is  the  turbulence  limited  coherence  diameter,  D is  the  tele- 
scope aperture  diameter*,  and  r0L  is  the  diffraction  limited  modulation 
transfer  function  associated  with  the  aperture.  For  the  case  where  h' s» 

If  then  since  |»f#  - k|  can  be  very  small  we  must  write 

<|t(h'-h)|2>t  ~ |tol(h'-  If) 1 2 expC-6.88  [|tf- S|/(k  r0)]*/»]  , 

for  Kf  as  K , (75) 


where  k = 2n/\  . To  accommodate  both  the  case  where  »?#  « h , and 

the  case  where  they  are  significantly  unequal,  we  write,  combining  Eq.  ’s 
(74)  and  (75) 

<|T(S'-!f)|8>T  « |TDL(^-H)|2exp(-6.88[|If'-7|/(kr0)]^3) 

+ 0.435  (r0/D)a  T#L(It/ - h)  . (76) 

We  argue  that  T(If)  is  approximately  a zero  mean  gqussian  ran- 
dom variable  so  that  odd  moments  will  vanish.  Accordingly,  we  may 
write 


* We  assume  here  a circular  clear  aperture. 


* 
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<T*(*'  + K)  T(h')  t(h)>t  *,  0 


<T*0?#  - *)  T(*')  T*(*)>T  » 0 . 


<t(k'  - H)  T*(h')  t(h)>t  0 , 


<t(k'  + *)  t*(JT/)  t*(*)>t  * o . 


(77a) 

(77b) 

(77c) 

(77d) 


Finally  we  argue  that  since  t(h)  is  a zero  mean  gaussian  random 
variable  the  fourth  moment  of  its  magnitude  should  be  equal  to  twice  the 
square  of  the  second  moment,  i.  e.  , 


<iT0?>r>,  = 2[<iT(U)|»>,r  . 


(78) 


Thus,  if  It*  and  h are  nearly  equal,  so  that  t(>?/)  and  t(^)  are  well 
correlated,  we  expect 


<|T(*')|S|T(*)|a>T  - 2<|T(>?')|a>T  <|T(*)|3>T  . 


for  H fts  H 


(79) 


On  the  other  hand  for  5?'  and  h significantly  different  in  value,  so  that 
t(>?* ) and  t(>?)  are  effectively  uncorrelated  we  get 


<|t(*')|3|t(*)|3>t  - < | t(h* ) |a >T  < |t(77) | 3 >t  . 


for  S'  * 5? 


(80) 


Based  on  the  information  we  have10 on  the  way  the  correlation  between 
r(S')  and  t(h)  decreases  with  increasing  value  of  |h'-  >t | , we  argue 


R 


I 
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that  a reasonable  form,  interpolating  between  the  results  of  Eq.  's  (79) 
and  (80)  is 

<|T(*')l3|T(S)|a>T  * <|T(K/)|a>r  <|t(h)|2>t 

X (l  +exp{-6.88[p-H|/(kr0)]^3))  . (81) 

Making  use  of  Eq.  's  (74),  (76),  (77),  and  (81)  we  can  now  re- 
write Eq.  (73)  as 

CL(*,*')  = N4  |l0(>?')|3  |Io0?)|8  (0.435)3(r0/D)«  tDl(x')  tol  (S) 

X exp{-6.88  [(S'  - >?|/(k  rQ)]e/3} 

+ N3  |l0(^+J?)|a  (0.435)  (r0/D)2  tdl  (*'+£) 

+ W 1 10  (k" - k) | s ff(0.435)  (r0/D)»Toi(H'-  *) 

+ | tol  (h/  - *)  | 3 exp{-6.88[|H/-  n|/(k  r0)]6/33  J . (82) 

With  this  result  in  hand  for  the  covariance  of  the  spectral  strength 
obtained  with  the  Labeyrie  method  of  spectral  interferometry,  we  can  cal- 
cluate  the  expected  error  in  this  process.  To  make  this  quantity  meaning- 
ful we  also  need  an  expression  for  the  expected  signal  strength.  We  can 

obtain  this  from  Eq.  (64),  which  in  conjunction  with  Eq.  (74)  leads  to  the 
result  that 

<spt,if))  = N2  |i0P)p  (0.435)  (r0/D)a  T0lP?)  . (83) 

In  the  next  section  we  use  Eq.  (82)  and  (83)  to  form  an  estimate  of  the 
accuracy  of  speckle  interferometry  measurements. 
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3.  9 Speckle  Interferometry  Measurement  Accuracy 


In  speckle  interferomotry  wo  form  a set  of  M statistically  in- 
dependent speckle  patterns,  and  then  by  fourier  transforming  and  aver- 
aging over  these  M patterns,  we  develop  an  estimate  of  S(J?,h).  The 
expected  value  of  this  quantity  is  (S(h,h))  , just  as  given  by  Eq.  (83)  . 
The  covariance  associated  with  this  quantity  will  be  Cl(h,m')/M  , with 
the  value  of  Cl(k,»?/)  given  by  Eq.  (82).  Nominally,  this  establishes 
a "signal-to-noisc  ratio"  in  the  spatial  frequency  domain,  SNRH  > 
which  we  can  write  as 

SNRk  = (S(h,  H))a/[CL(K,T?)/M]  . (84) 

Making  use  of  Eq.  ’s  (82)  and  (83)  with  appropriate  approximations,  we 
can  rewrite  Eq.  (84)  as 


M 

hNR>t  = 1+[N  |ro6)p(0.435)(r0/D)=  T#l(*)ra  * (85) 

where  we  have  assumed  that  r0  is  very  much  less  than  D . The  inter- 
pretation of  this  formula  is  in  a sense  straightforward.  If  the  average 
photon  detection  events  per  single  short  exposure,  N , is  large  enough 
so  that  N 1 10(~)|  s(0.  435)  (r0/D)  tdl  (h)  is  significantly  greater  than  unity, 
then  the  signal-to-noisc  ratio  will  not  be  improved  by  any  increase  in  the 
photon  flux  — the  process  will  not  be  limited  by  shot  noise  but  rather  will 
be  entirely  limited  by  turbulence  effects  and  the  need  to  average  over  many 
realizations  of  the  random  turbulence  effects. 

Unfortunately,  this  spatial  frequency  domain  signal-to-noisc  ratio, 
while  it  is  useful  in  providing  us  with  some  insight  into  the  required  photon 
detection  event  rate  is  not  directly  relatable  to  the  accuracy  of  a measure- 
ment quantity  of  interest.  It  is  measurements  in  the  spatial  (rather  than 
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the  spatial  frequency)  domain  that  are  basically  of  interest  to  us.  It  is 
in  terms  of  the  evaluation  of  this  type  of  quantity  that  we  wish  to  deter- 
mine the  measurement  accuracy  and  the  associated  signal-to-noise  ratio. 
A s an  example  of  this  we  consider  the  observation  of  a binary  star  pair. 
We  consider  a binary  star  pair  with  an  angular  separation  © and  with 
an  intensity  ratio  p (with  p less  than  unity).  As  an  example  to  be 
pursued  here,  we  ask  how  accurately  we  can  determine  the  value  of  p 
with  a set  of  M short  exposure  photographs  if  there  are  an  average  of 
li  (l  + p)  detected  photon  events  per  short  exposure. 

For  this  pair  of  binary  stars,  making  use  of  Eq.  's  (2)  and  (3)  we 
can  see  that 

I0(r)  = (1  + p)"1  [6(?)  + p 6 C?  - 0)]  . (86) 

so  that 

10(k)  = (1  + p)"1  [1  + p exp(-  i k • 1)]  . (87) 

A set  of  M short  exposures  will  be  used  to  form  the  quantity  S(h)  which 
estimates  (S(k,h))  . For  observation  of  turbulence  statistics,  based 

nost  likely  on  images  of  a bright  single  star,  a normalization  factor 
sn(h)  will  be  formed,  which  will  have  a nominal  value  of 

Mh)  = 0.435  (r0/D)»  tol(h)  , (88) 

for  large  values  of  |h|  , but  actually  corresponding  everywhere  to  the 
basic  modulation  transfer  function  of  the  speckle  interferometry  process. 

The  basic  speckle  interferometry  process  would  proceed  by  form- 
ing the  ratio  5(T?)/sm(h)  , which  in  accordance  with  Eq.  183)  and  (87) 
should  approximate  the  power  spectrum  of  the  binary  star  image,  i.  e.  , 
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Making  use  of  Eq.  (87)  this  takes  the  form 


S(h)/sn(x)  m Na(l  + Prs[l  + pa  + 2p  cos(#  • ?)]  . (90) 

This  quantity  would  then  be  fourier  transformed  to  form  an  estimate  of 
the  correlation  function  associated  with  the  binary  star  pair's  image. 

We  obtain  for  the  correlation  function 

J (?)  = (2tt)”2  Jd>?  exp(i  h • 7)  "S(h)/sn  (£)  . (91) 

If  we  substitute  Eq.  (90)  into  Eq.  (91)  and  carry  out  the  fourier  trans- 
formation integration  we  get  as  the  nominal  correlation  function 

J(f)  = N2(l  + p)”2  [(1  + p2)  6(r)  + p J(r  - ?)  + p6(r  +0)]  , (92) 

which  is  exactly  the  correlation  function  we  would  expect  for  the  binary 
star  pair's  image  — a main  lobe  of  relative  strength  1 + p2  at  the  origin, 
and  two  side  lobes,  each  of  relative  strength  p at  displacements  of  ±?  . 
In  writing  Eq.  (92)  we  have  used  the  equivalence 

0(r)  = (2n)“a  /dr  exp(i  >?  •”?)  . (93) 

Since  the  main  lobe  and  side  lobes  can  be  separately  measured  we  can 
form  the  ratio  R whose  nominal  value  is 

- P/d  + Pa)  • <94) 
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but  whose  measured  value  is  R = R + bR  . This  allows  us  to  calculate  the 
relative  intensity  of  the  two  binary  stars. 

Solving  for  p we  get 


1 -HI  - 4/fa)1/a 
9 = 2R 


Qur  basic  concern  is  with  the  effect  of  some  fractional  error 

« = bR/R  , <9 

in  the  actual  determination  of  R (due  to  the  statistical  nature  of  the 
speckle  interferometry  process)  on  the  fractional  error 


p,  = 6 p / p . *97' 

in  the  determination  of  the  relative  intensity,  p , of  the  pair  of  binary 
stars.  To  obtain  this  relationship  we  write,  in  accordance  with  Eq.  (95), 


p + 6p  = 


\3  T 1/3 


2(*  + bR) 


This  can  be  rewritten  as 


p + 6p  = 


i + n - 4*s  (i  + «)8i^s 


2/?(l  + «) 


which  we  can  approximate,  for  c small,  by  the  expression 
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then  the  value  of  % is  obtained  when  "5(h)  in  Eq.  (91)  is  replaced  by 
(S(H  ,h))  , and  that  substituted  into  Eq.  (92),  we  see  that  we  can  write 


6<p  = (2TTH>,rd*  expUa.?)[S(>?)  - <S(k.7)>]/sn(7) 

(2TT)-a  Jd*  <S(*  ,?))/sn(h) 

_ (2n)~a  fdg  exp(i  * . ?)  <S(k  , k)>/sn(k) 

(2tt)“s  Jdx  <S(k  , k))/sn(h) 

x (2nra  ,fc W [SQ?)  - (S(h  , %))]/ sN (k)  . (1Q4) 

(2n)“a  J*d>?  <S(h  , k))/sn(h) 

We  note  that  since  sN  (h)  is  a normalization  function  chosen  to  every- 
where correspond  to  the  speckle  process's  transfer  function,  then  in 
accordance  with  Eq.  (83)  , [and  taking  note  of  Eq.  (88)],  we  can  write 

<S(J?,S)>/sN(#)  = T7a;d*  1 10 (7?) 1 3 

= T7a  C . (105) 

It  is  easy  to  see  from  consideration  of  Eq.  (87)  that 

C = JcPt  |l0(Z)|a 

= ( 1 + p)“a  Jdx  | 1 + p exp  (-  i h • "8)  |a 

[O  + P3)/(l  + P)a]  n (kD)a  , (106) 

where  the  factor  of  tt  (kD)a  appears  here  as  a consequence  of  the  fact 
that  the  range  of  the  H-integration  is  a circle  within  which  |h  |<  kD,  i.  e. , 
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less  than  the  diffraction  limited  transfer  function  cut-off  spatial  frequency, 
kD?  Noting  that 

o - (2n)~3  -fd*  exp(i  h •?)  <S(?T,  h))/s»,(h) 

' „ _ # 1 1 07) 

(2nra  /<fi?  <S(S,Z)>/sN(S) 

we  can  now  write  in  place  of  Eq.  (104) 

. - jst  [s*p» ^ [30?)  . <spf  .*)>]  . U08) 

Our  interest  is  in  the  evaluation  of  the  mean  square  value  of  the  fractional 
error  in  c , the  measurement  of  the  ratio  of  the  two  peaks  in  the  correlation 
function. 


°l  - <l«|-> 

= (hpQf r®  f fd*  dH'  r ex£(-  exp(iH  .?).Ar 

L JL  s„(>?)  J 

x <[S(h')  - <S (h',h'))]*[5(h)  - <S<r  ,k)>3>  . (109) 


In  as  much  as  S(h)  is  the  measurement  approximation  to  S(h,h)  except 
averaged  over  M short  exposures,  then  we  can  replace  the  ensemble  aver- 
age in  Eq.  (109)  by  C^k.kO/M  , where  Ct  is  defined  initially  in  Eq.  (65), 
and  in  the  form  we  shall  use  in  Eq.  (82).  This  allows  us  to  write 


= (n ac*r»M-i ;/  d^  ^ 


x"{W1rf.^|-8|fD(1f)|B(0.435).(r0/D)*  TDt(*')  tol(k) 


* Assuming  0 > \/D  we  have  dropped  a term  of  the  order  of  pJ  (2n6D/X)/ 
(2tt0D/\)  in  Eq.  (106)  . 
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X expt-6.88  t|x'-HlMkre>r»] 

+ T7i. * 3  |I0(h#+  K)JS(0.435)(rQ/D)a  V<a+*) 

+ TTa  |10(h  - h')|3  ^0.435)(r0/D)2  t01(k#  - *) 

+ 1 tpli*  - *')\2  eXP  I-6*88  U*  " *'|/<kro>36/3}  ^ ^ ' 

Making  use  o i Eq.  (88)  we  rewrite  this  in  the  form 

x ll0(x)|a  exp[-6.88[|x'-  n|/(kr0)]5/  } 

[Ipf^  +H*)13  T.,(h  

+ T73(0.435)(ro/D)2  t0L(x)  tol(h') 


^ lln(H  -^)la  T 

+ "H3  (0.  435)(r0/D)3  tol(k)  t6L(x#) 


+ U~.(>?  -^)lalT0L(x  -x')l3exp{-fe.88[)H  -gM/^r0fl^j 

+ T7a(0.435)8(ro/D)«  Tet(H)  tol(h')  ^ 


In  these  integrations,  though  it 
stood  that  the  integration  range 
the  diffraction  limited  transfer 


i.  not  explicitly  stated,  it  is  to  be  under- 

only  include,  value.  of  » and  * for  which 

function,  V • doe‘  "0t  Va““h' 
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Rather  than  attempt  to  exactly  evaluate  the  multiple  integral  in 
Eq.  (Ill)  it  is  convenient  to  approximate  the  result  in  the  two  limiting 
cases  of  1)  very  large  photon  flux,  i.  e.  , N very  large,  and  2)  a very 
weak  photon  flux,  i.  e.  , TT  very  small.  For  the  first  case  we  can  write 

= o,a  . for  'R » 1 

(WMH  d *'  [exp(-i  *'•*)  - **] [exp(i  Z . ?)-*>] 

^ ly^)|a  |I0  W|Sexp£-6.  88L|k'-  h |/(k  r0)]»/»)  . (112) 

If  we  now  make  use  of  Eq.  (87)  and  furthermore,  restrict  our  attention 
to  the  case  of  binary  star  pairs  whose  angular  separation,  1*1  , is  too 
small  to  be  resolved  by  conventional  imagery,  i.  e.  , 

|?|  « (kr0)“l  , (113) 

then  we  can  make  the  approximation  that 

°«s,i  (1  -«■>) 

X J*dvl  exp{-6.  88  [|il  |/(k  r^]6^3}  . (114) 

It  can  be  shown  that 


J*du  exp  £-6.88[|Tl|/(kr0)]®/3  3 = 2n  C(6.  SS)”3'3  k r#  30 


= 0.3726  (k  r0)a  T (6/5) 


= 0.3422  (krQ)2  . (115) 

Now  making  use  of  Eq.  's  (95)  and  (106)  to  provide  values  for  R and  C , 
respectively,  and  using  Eq.  (115),  we  can  reduce  Eq.  (114)  to  the  form 


o 


2 

e,l 


0. 1089 
M 


(r0/D)2  (3± 


(116) 


Making  use  of  Eq.  's  (102)  and  (95)  we  can  see  that  in  this  case  the  rms 
error  in  the  measurement  of  the  ratio  will  be 

ffp.i  = P2  <“i> 

= P2*u2x 

pS  (1-4^  F2 

= p2  (t~T^)(2P  ~ ^ ^)(Q-  1089)(r0/D)g/M 
= (2p  - p2  + 2p2)  (yt4)  (0.1089)(ro/D)2/M  . (117) 

This  represents  our  basic  result  for  our  ability  to  measure  the  relative 
intensity  of  a pair  of  (unequal  intensity)  binary  stars  for  case  1,  when  the 
photon  flux  rate  is  very  high  and  the  error  is  associated  with  the  imper- 
fect average  over  turbulence  effects  provided  by  M short  exposures. 

- 72  - 


For  case  2,  where  the  photon  flux  rate  is  very  low  we  can  approx- 
imate a€3  as  given  in  Eq.  (Ill)  by  the  expression 

ai,a  RS(C2/?2M)"^*J'dx  d*#  [exp(-i  H,*!f)-^’*][exp(i  h • 

x |l0(H  + y>|3  TDl(g  + >?')  + |I0(H  - 3')|2  tdl(>?  - 
N3  (0.435)  (r0/D)3  T0L(*)  TOL(>t0 

m 2( C3* 8 M)"i  JJdn  d*' R*  ll0(H-y)|  tdl(h  - kQ  . 

TTs(0.435)(ro/D)3  tolO?)  tol(S') 

(118) 

Making  use  of  Eq.  (87)  and  approximating  the  diffraction  limited  transfer 
function  r0L  , (whose  value  varies  from  1 to  0)  by  the  constant  § , we  can 
approximate 

o3a  « 2( C8 M)"1  L(  1 + Pa)/il  + p). J $ [n(kD)a^8  . (H9) 

r73(0.435)(r0/D)8(i)8 

Now  making  use  of  Eq.  's  (95)  and  (106)  we  can  reduce  this  expression 
to  the  form 


a 13.8  /D\3  (1  + p)3 

ae,3  ~ MW3"*  ( r0  ) ( 1 + p2) 


(120) 


Following  the  procedure  of  Eq.  (117)  we  see  that  in  this  case  the  rms  error 
in  our  measurement  of  the  ratio  of  the  intensity  of  the  two  stars  in  the  bi- 
nary star  pair  will  be 


o 


3 

p,  3 


a 

3 0-0 

p II  S — 

(1-4/f8)!/3 
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(121) 


= oa 


1 + pa  1 

- 13.8  / 

D Na 

(1  + o)a  > 

1 - pa  1 

.MN5  V 

TJ 

(I  + Pa)  r 

8 

f JLLfi. 

a 

N®  V r0  . 

) 1 - p 

P 

• 

It  is  interesting  to  note  here  that  since  W is  proportional  to  Ds , in- 
creasing the  aperture  area  will  reduce  the  rms  error  in  our  measured 
value  of  p , but  only  as  the  first  power  of  the  aperture  diameter. 

It  is  convenient  to  combine  our  results  for  the  two  limiting  cases 
of  a very  large  photon  flux  rate  and  of  a very  low  photon  flux  rate, 
writing 


3 - 


CS  + a 3 
p.i  + cp.a 


= s4r2(^)s 


+ 8 r D \s  1 + p 

' r J 1 - p 


M7Ta  v 


(122) 


With  this  result  established  we  are  now  ready  to  consider  the  magnitude 
of  noise  effects  in  the  Knox-Thompson  speckle  imagery  technique.  We 
take  this  up  in  the  following  sections. 


l 


. 


3.  1 0 Accuracy  of  Speckle  Imagery 


For  consideration  of  noise  effects  in  speckle  imagery  our  basic 
concern  is  with  error  in  the  estimation  of  (Stf,*'))  for  >?'  approximately 
but  not  exactly  equal  to  h . In  the  speckle  imagery  process  the  phase  of 
<S(it  .x7))  is  determined  for  an  array  of  values  of  >t  and  Tt*  that  "cover" 
the  7?-space  range  of  interest.  These  phase  values  may  be  considered  to 
be  equal  to  the  difference  in  the  phase  shift  associated  with  the  image  pattern's 
spatial  frequency  components  IQ  (h)  and  1 0 ( 7?' ) . By  a least  square  method 
these  phase  differences  are  then  combined  so  as  to  span  the  entire  spatial 
frequency  domain,  yielding  the  (absolute)  phase  shift  to  be  associated  with 
each  spatial  frequency  component  in  image  patterns  spectrum,  10(h)  . It 
can  be  shownnthat  if  the  mean  square  error  in  the  determination  of  the 
phase  difference  is  constant  over  the  spatial  frequency  domain,  and  is 
equal  to  cAC p , then  the  mean  square  error  in  the  calculated  phase  will  be 

cr*  * f [1  + {»)]  oA%  , (123) 

where  91  is  the  number  of  points  in  the  spatial  frequency  domain  on  which 
the  phase  is  to  be  determined,  and  over  adjacent  pairs  of  which  the  phase 
difference  is  measured.  Knowing  the  values  of  a ^ , the  mean  square 
phase  difference,  and  of  51  , the  number  of  measurement  points  in  the 
spatial  frequency  domains,  we  can  then  consider  the  mean  square  error 
in  our  knowledge  of  the  phase  of  each  of  image  spatial  frequency  components, 
o^p  , to  be  known.  We  may  take  as  a measure  of  the  quality  of  the 
speckle  imagery  ultimately  recovered,  nominally  estimating  the  effective 
resolution  as 


(124) 
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*%\  = exP  (-  O*)  . 


This  may  be  considered  to  be  a measure  of  the  Strehl  definition  of  the 
recovered  image  relative  to  the  ideal  diffraction  limited  image. 

Our  immediate  concern  is  with  the  evaluation  of  . We  take 

this  up  in  Section  3.  11,  recognizing  that  the  measured  phase  difference 
is  actually  a function  of  k , and  so  busy  ourselves  there  with  calcu- 
lating 0^(k)  . In  Section  3.  12  we  introduce  consideration  of  the  entire 
spatial  frequency  domain,  and  evaluating  31  , so  that  we  can  then  deter- 
mine ct„  and  ^S|  . We  shall  develop  both  of  these  quantities  in  Section  3.  12. 
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3.11  Speckle  Imagery  Phase  Difference  Error 


The  starting  point  for  our  calculations  here  is  our  ability  to  evaluate 
( |S(>?  ,h7)|  s)  . From  this  we  can  evaluate  the  mean  square  error  contained 
in  our  estimate  of  (S(h  , h7))  as  formed  from  M frames.  Since  S(h  ,k7)  is 
a complex  number  the  error  is  as  likely  to  be  in  the  real  part  as  in  the 
imaginary  part  — or  put  in  more  relevant  terms,  as  likely  to  be  parallel 
to  (S(#,k7))  as  it  is  to  be  perpendicular  to  (S(5t  ,h7))  in  the  complex  plane 
containing  (S(h,k7))  . Thus,  based  on  the  assumption  that  after  averaging 
M frames  the  error  will  be  small,  we  can  approximate  the  mean  square 
error  in  our  estimate  of  the  argument  of  (S(h  , 5?7))  , which  is  the  mean 
square  phase  error,  as 

.V0t>  - M-1  l<gff.*»|‘3  . (U5) 

l<sp:  ,;')>|a 

Making  use  of  Eq.  (59)  we  can  write 
<|S(S  ,S')|S>  = <|R*(h')  R(h)-R(h  - 5?7)|8> 

= <R*(>?/)  R(Jf7)  R*(h)  R(h)> 

- <R*(h7)  R(x)  R*(h  - h')) 

- <R*(tt)  R(h')  R(k  - It7)) 

+ (R*(5?  - %')  R(Z  - h7))  . (126) 

■ , < jolting  the  formalizm  developed  in  Section  4 and  substituting  as 

,n  . using  Eq.  ’a  (37)  and  (38a,  b,  c,  and  d),  we  can  rewrite 

I tit)  •• 
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<|S(h,h')|3>  = <M4(-  h'  , - H , h)>t 

- <M3(-S',  * .H'-H)),-  <M 3(-»?  -S')>T 

+ >?,  7t  - *')>T  • (127) 

Now  substituting  Eq.  's  (45),  (50),  and  (52)  into  Eq.  (127)  we  can  obtain 
the  result  that 

<|S(>t  .tf)!3)  = (TT4<I(-ff)  !("')  I(-H)  I(5?)>t 

+ TT3<I(0)  f(-5?)  1(h)>t  + N3<I(-S'->?)  l(K')  1(Z))T 
+ N3  <1(k  - >?')  1(5?')  I ( -h))t  +TT3<I(h/-  k)  I(-h')  I (7?)>t 
+ T73  <I(h'  + h)  I(-*')  1(-H)\  + N3<I(0)  I(-h')  I(h')>t 
+ N3  (1(0)  I(0)>T  + N3  <I(-h'-h)  I (k7  + k)>t 
+ N3  <l(Jf  - h')  I(*'-  h)>t  + N8<I(-h)  I(h)>t 
+ Tia  <i(h)  i ( ->?)>T  +TT3(I(-^)  i(h')>t 

+ TT3<i(»?')  f(-*')>T  + TT  <I(0)>T  3 

- (TT3<I(-h/)  1(h)  I(k'-  *)>t  + NS<I(h  - h')  for-  h)>t 
+ Na  <I(-m)  I(>?)>t  + NS<I(*')  1(-h/)>t  +N(1(0)>t) 

- (N3  <!(->?)  T(h')  I(*  - h')>t  + N3<1(h/- It)  I(h  - h#)>t 
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+ TJa  <!(-*')  I(H/)>r  + Na<I(Tf)  !(-*)>,  + N <l(0)>r  ) 


t 


+ {Na<i'(*'- h)  1(3  - *')>T  + N <1(0)>T  ] . (128) 

Cancelling  terms  as  appropriate,  this  expression  can  be  reduced  to  the 
form 

<|S(H  , = N4  <1(-^)  1(3')  !(-*)  I(7?)>t 

+ 7I3<1(0)  I(-k)  1(3)>t  + N3<l(-K#-i?)  I(3')f(3)>T 

+ TT3  <!(*'  + 3)  r(-K;)  i(-k)>t  + n3<i(0)  M-S')  f(3')>T 

+ N8<1(0)  I(0))t  + TJa  <!(->?'-  k)  1(h' + h)>t  . (129) 

Now  making  use  of  Eq.  (7)  and  removing  from  inside  the  angle  brackets 
all  the  factors  which  do  not  depend  in  a statistical  way  on  the  turbulence 
induced  wavefront  distortion,  as  well  as  using  Eq.  's  (56)  and  (57),  we 
can  rewrite  Eq.  (129)  as 

( | S (x  , 3' ) | a)  = N4|f0(S)|3  |10(*')|8  <|t(3)|3  |r(3')|a> 

+ N3CIo(0)  |r0(7?)|c  <T(0)  |r(TT)r>T  + I*(3  +H')l0(^)l0l7)<T*(Tf+^)T(»f/)T(Jf)X 

+ r0(>?  + ^)  r0V)  i0*(*?xt(^  +r)T*(H0T*(*)>t 

+ io(o)|i(«0|2<t(o)|t(T?/)|8>t) 

+ TTat|io(0)|°<|T(0)l3)T  + |f0(3  + 3')|»<|t(3 +*')|»>T]  . (130) 


i- 

I 


I 1 


I 
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Taking  advantage  of  the  fact  that  Io(0)  and  t(0)  are  each  equal  to  unity, 
this  can  be  further  reduced  to  yield  the  result  that 


<|s(h.*')|3>  = N4|r0(S)|»|i0<*)l8  <|T(*)|a|T(S/)|a> 

+ ^3Uf0(^)|a<|T(«)|a>T  +i*(5?+^r0(H)i0(>to<T*a?  +^>t(«)t(^)>t 

+ I0(*  + *')f?(*)rj<*')<T(*  + kO  t*(>?)  t*(h/))t 

+ |Io(*')|a<|T(3')|a>T) 

+ TTa[l  + |I0(W  + *')|8<|T(*+S')|a>T3  • (131) 

From  Eq.  (64)  we  see  that  we  can  write 

I <s<*,S')>|s  = N4 |r0(J?)la|r0(^)la|<T*(H')T(H)>|3 . (i32) 

Making  use  of  Eq.  (74)  we  can  rewrite  this  as 

l<S<*  • *)>P  = TT*|I0(^)|3  |ro(^)|a(0.435)3(r0/D)*T(>L(7?)TDL(>?0  . 

(133) 


Similarly  making  use  of  Eq.  's  (74),  (77),  and  (81)  along  with  Eq.  (131) 
we  can  write 

<|S(K  ,*')|a>-|<S(*  ,*')>|a  = N4  | ro(n)  | 3 |ro(>t/ ) | a ( 0.  435)3  (r0/D)4 
* t6l(k)  Tou(ff')  exp(-6.  8 8 C l ~ *|/(kr0)]6'33 

+ N3{|lo0?)|a(0.435)(ro/D)=TDL(Tf)  + |l0(*0|a(0.435)(r0/D)*TDt(S')) 


+ Na[l  + |I0(>?  + *')|3(0.435)  (r0/D)*TDL(K  + *')}  . (134) 


Substituting  Eq.'s  (133)  and  (134)  into  Eq.  (125)  and  agreeing  to  restrict 
attention  to  values  of  h and  k '*  sufficiently  close  together  that  we  may 
approximate 

tdl(h)  « tol (h;)  , (135) 

and  for  a nominally  smooth  image  spectrum 

r0(*)  * I0(*')  , (136) 


we  get 

aA%(H)  = (2M)n  2exp{  -6.  88[  [h7-  h |/(kr0)]6/3  ] 

+ 2[N|IQ0?)|2(0.435)(ro/D)2  Tot(jt)]-x 
+ [N|I0(^)|2(0.435)  (r0/D)»  tdl(h)]“= 

+ I0(h  + S')  rDL(^  +T?')[TT|l0(H)|4(0.435)(r0/D)2|ToL(T?)|2]-i2 

(137) 

Based  on  the  assumption  that  h and  h * are  close  enough  that 

6.88[|>?  - *'|/(kr0)]6/3  « 1 . (138) 

it  is  reasonable  to  approximate  Eq.  (137)  as 


- 81  - 


o*  (*)  - (2Mp[l  + W|l0(H)|3(0.435)(r0/D)aTBL(»?)ri}3  . (139) 


The  detected  photon  flux  requirement  for  this  result  not  to  be  photon  shot- 
noise  limited  is  that 

TT  » [|lo0?)|a  (0.435)  (r0/D)3TDL(S)]  . (140) 

Clearly,  as  we  direct  our  attention  to  higher  and  higher  spatial  frequencies 
where  tdl(h)  is  smaller  and  smaller,  and  where  presumably  so  is  the 
image  spectrum  we  require  a higher  and  higher  detected  photon  flux  for 
our  results  to  not  be  photon  shot-noise  limited.  No  matter  what  the 
photon  flux  rate,  at  some  spatial  frequency  performance  will  become 
shot  noise  limited. 

With  Eq.  (139)  in  hand  we  may  now  turn  to  our  ultimate  task,  the 
evaluation  of  the  resolution  associated  with  the  Knox-Thompson  speckle 
imagery  technique  for  a finite  set  of  short  exposures.  We  take  this 
up  in  the  next  section. 


\ 
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3.  12  Speckle  Imagery  Resolution 


The  basic  equations  for  evaluation  of  the  achievable  resolution  in 
speckle  imagery  are  Eq.  's  (123),  (124),  and  (140)  . In  addition  to  these, 
we  need  expressions  for  91  and  some  indication  of  how  f0(x)  and 
tdl  (*)  may  be  expected  to  vary  as  a function  of  >?  . We  may  approximate 
the  form  of  tol  (x)  by  the  expression 

t0l(2)  « 1-1.25  [|H|/(kD)]  + 0.25  C|K|/(kD)]4  , (141) 

where  k = 2tt/X  and  D is  the  telescope  aperture  diameter.  As  a prac- 
tical matter,  over  most  of  the  spatial  frequency  range  of  interest  we  may 
simply  write 

tol(h)  w 0.  5 , (142) 

as  a sufficiently  good  approximation  for  our  purposes. 

If  the  object  being  viewed  has  an  angular  extent  ©fl  , and  is  rea- 
sonably "isotropic"  in  its  size  and  complexity,  and  is  complex  in  appear- 
ance (i.  e.  , contains  many  features  of  all  sizes),  then  it  is  probably  rea- 
sonable to  approximate 

|io(*)|a  = Cl  +(|H|/H0)a]-1  , (143) 


where 


*o  = n 


(144) 


For  the  spatial  frequencies  of  interest  we  can  approximate 


- 83  - 


x k i i r 1*11  if  rTiirrifciai^- — — 


im*>i2  « k/|h|)2 


(145) 


To  estimate  Dt  it  is  convenient  to  introduce  the  factor  £ relating 
the  spatial  frequency  increment,  i.  e.  , |»f  - h'|,  in  our  speckle  imagery 
process  to  krQ  . If  the  object  size,  ©0  , is  small  enough  so  that  the 
turbulence  limit  sets  the  apparent  image  size  and  the  extent  of  each  of 
the  recorded  speckle  patterns,  then  we  know  from  analysis  of  the  Knox- 
Thompson  algorithm  that  we  must  use  a spatial  frequency  separation 
such  that 

|k-  *'|  * p kr0  , (146) 

where  ^ is  a number  less  than  one-half  but  probably  greater  than  one- 
tenth.  Nominally  we  take  ^ = 0.25  for  this  analysis.  If,  however,  ©0 
is  large  enough  so  that  the  objects  actual  size  sets  the  extent  of  the  in- 
dividual speckle  patterns,  so  that  it  is  v*1  times  as  large  as  it  would 
be  for  a purely  turbulence  limited  extent,  (where  v is  a number  less  than 
unity),  then  the  lowest  spatial  frequencies  of  interest  would  be  of  the  order 
of 


*0  = "/®o 


= i v k r0 

Thus  we  may  approximate 

{i  v , for  0O  > \/r0 

V-  , for  0Q  < \/r0 


(147) 


(148) 
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and  write 


In  estimating  the  image  phase  error  at  spatial  frequency  k we 
shall  make  the  convenient  assumption  that  we  are  not  going  to  bother  to 
obtain  data  for  spatial  frequencies  whose  magnitude  is  greater  than  |3| 
Accordingly,  there  will  be 


31  = TT  L|*|/(kr0)] 


It  is  convenient  to  note  that  in  general  we  may  consider 
significant  fraction  of  the  optics  spatial  frequency  cutoff 


where  o is  a number  nominally  between  one-tenth  and  unity.  Accord 
ingly,  we  may  rewrite  Eq.  (150)  as 


Thus 


In  general,  we  may  expect  the  quantity  in  the  square  brackets  to  be  domi 
nated  by  the  D/r0  dependent  term,  so  that  we  can  approximate 


It  is  obvious  that  the  D/rQ  dependence,  while  not  actually  negligible 
will  have  only  a minor  impact  on  the  overall  phase  error,  a * . Doubling 
the  telescope  aperture  diameter  may  effect  o£  through  its  effect  on 
°<f<p  . but  will  have  only  a very  slight  effect  on  through  the  fa(D/ rQ)- 
term. 

To  evaluate  we  substitute  Eq.  (145)  into  Eq.  (139).  This 

allows  us  to  write 

»*%<*>  * (2MH  {1  + [0.435  N (r0/D)8(H0/|S|)»  tol(h)]-1]8.  (155) 

Making  use  of  Eq.  's  (142),  (144),  and  (151)  we  can  reduce  this  expression 
to  the  form 

»*%(*>  « (2Mp|l  + {0.435  N(ro/D)3[(n/0o)/(akD)]3(i))“i|.S 

- (2Mrl{l+tO-05«-*(r0/D)=[(i^-)SNj'‘]-»|3  (J56) 

The  quantity  [(X/D)/@]s  can  be  considered  to  be  the  number  of  resolution 
elements  in  the  target  object  (whether  or  not  the  target  is  greater  than  or 
smaller  than  \/rQ)  , and  so  the  quantity  in  the  square  brackets  in  Eq.  (156) 
can  be  considered  to  be  the  number  of  detected  photon  events  per  exposure 
time  per  resolution  element.  It  is  this  rather  than  any  other  measure  of 
the  total  target  brightness  that  determines  whether  or  not  the  speckle 
imagery  process  is  or  is  not  limited  by  shot  noise  effects.  The  basic 
criteria  for  the  process  not  to  be  shot  noise  limited  is  that 

N (3^)  > 20«a (D/r0)2  • (157) 
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Since  TT  will  only  increase  as  Ds,  it  follows  from  Eq.  (155)  that  the 
larger  the  aperture  diameter  the  brighter  the  target  object  has  to  be  to 
yield  a speckle  imagery  result  that  is  not  limited  by  photon  shot  noise. 

If  the  target  object  is  sufficiently  bright  that  aACp  is  essentially 
independent  of  N , then  the  mean  square  phase  error  in  the  recon- 
structed target  image  using  speckle  imagery  techniques  is,  in  accord- 
ance with  Eq.  's  (154)  and  (156)  , 


* . fr(D/r0) 

6 M 


(158) 


and  the  associated  resolution  is,  in  accordance  with  Eq.  (124) 


= (r0/D)a/s«)  . 


(159) 


When  the  target  brightness  is  small  enough  that  photon  shot  noise  effects 
are  significant,  then  we  can  write 


a 


3 

<P 


J 


&i(D/r0) 


(160) 


where 


Y 


I N/£tt®^ 

L iirD* 


J [“  (f  ) (rj  j 


3 “1 


(151) 


We  note  that  Y is  a measure  of  the  target  brightness,  (corresponding  to 
the  first  square  bracket  quantity)  relative  to  the  brightness  required  so 
that  the  speckle  imagery  process  of  the  spatial  frequency  |k|  = ak  D 

will  not  be  limited  by  shot  noise  effects.  It  is  interesting  to  note  that 
this  required  brightness  level  is  independent  of  aperture  diameter,  D , 
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if  we  hold  the  spatial  frequency,  |k  J , constant,  but  increases  with  D3 
if  we  ask  that  the  spatial  frequency  always  be  a constant  fraction  of  the 
aperture  cutoff  frequency,  kD  — i.e. , if  we  hold  a constant  as  we 
vary  D. 

Substituting  Eq.  (154)  into  Eq.  (124)  we  obtain  as  our  final  result 
the  expression 

= (VDJCx+Y’V/tes)  , (162) 

the  speckle  imagery  resolution.  This  represents  the  last  of  the  results 
we  undertook  to  develop  in  this  chapter. 
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